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Two years ago it was pointed out! that the electrostatic valence rule 
requires that the aluminum ions in alkali aluminum silicates have the 
coordination number 4 rather than 6, oxygen ions in the crystal then 
being common to a silicon tetrahedron, an aluminum tetrahedron, and 
one or more alkali-ion polyhedra. Since that time examples of such struc- 
tures have been found. Thus the cubic crystal sodalite,? Na,Al;Si;0,.Cl, 
has a framework composed of tetrahedra, each of which shares its corners 
with four others, silicon and aluminum tetrahedra alternating. The 
three parameters determining the position of the oxygen atoms were 
evaluated from the x-ray data, and it was found that both the silicon 
and the aluminum tetrahedra are approximately regular, with the dis- 
tances Si-O = 1.60 A and Al-O = 1.74A. The sodium and chlorine ions 
are arranged in cavities in the tetrahedral framework. 

We have found that many other alkali and alkaline-earth aluminosilicates 
have similar tetrahedral structures. Three of these, shown by natrolite, 
the scapolites, and davynite-cancrinite, are described in this communica- 
tion. 

Natrolite—The zeolite natrolite, NazAlSi;0i).2H:O, forms ortho- 
rhombic pseudotetragonal crystals with axial ratios a:b:c = 0.97852:1:- 
0.35362. Data from oscillation and Laue photographs lead to the unit 


a=18.19A, 6 =1862A, c = 658A, 


containing SNaeAlSijs0i0.2H,O. The lattice is face-centered, and the 
space-group is shown to be C}} by the absence of certain prism plane 
reflections. 

The nature of the framework of the crystal can be deduced with some 
rigor from this information and the assumptions that aluminum and silicon 
tetrahedra of approximately the size of those found in sodalite exist in 
this crystal, and that they share corners but not edges or faces.* 8Si 
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FIGURE 1 


The natrolite string. Each 
tetrahedron has an oxygen ion 
at each corner and a silicon or 
aluminum ion at its center. 
Below is shown the projection 
of the string on (001). 
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must then be placed in Wyckoff’s position 8a, and the parameter can be 
taken as zero without loss of generality. Thus there are tetrahedra at 


000, 1/21/20, 1/20 1/s, 0'/2'/2, VW B/ s/s, 3/43 /a4/4, 3//8/4, 1/,3/,3/,. 


The distance from 000 to !/41/4!/4 is 6.67 
A, too far to be spanned by a single tetra- 
hedron. Hence two silicon and two alumi- 
num tetrahedra must be attached to each 
of these silicon tetrahedra. In order for 
each tetrahedron to share all its corners it 
is necessary for each of these groups to form 
a string with those obtained from it by 
the translational operation c, as shown in 
figure 1. The strings are rigid, and the 
dimensions given above for the tetrahedra 
lead to c = 6.57 A, in agreement with the 
observed 6.58 A. The strings combine by 
sharing tetrahedron corners as shown in 
figure 2. 

The structure as a whole is not rigid; it 
collapses, the strings rotating about their 
axes until certain minimum interionic dis- 
tances are reached. If all the tetrahedra 
were the same (silicon tetrahedra, say) the 
structure would have tetragonal symmetry. 
The introduction of the slightly larger 
aluminum tetrahedra reduces the symmetry 
to that of the orthorhombic space group C}}. 
With the angle 6 shown in figure 2 equal 
to 20° 50’, the assumed tetrahedron dimen- 
sions give calculated values a = 18.22 A, 
b = 18.56 A, in good agreement with the 
observed a = 18.19 A, b = 18.62 A. 

The sodium ions and water molecules 
occupy certain positions between the strings, 
such that the minimum Na-O and H,O-O 
distances are about 2.6 and 3.2 A, respec- 
tively. The channels along (001) are of 
such a size that on heating the crystal water 
molecules can escape along them without the 
rupture of the framework. In this way our 
determination of the structure of a zeolite 


leads to a simple explanation of their most characteristic property. The 
phenomenon of base-exchange is similarly accounted for; the sodium ions 
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can run along the channels and out of the crystal as other cations run in 
to replace them. 

The monoclinic pseudotetragonal zeolite scolecite, CaAlSis0,.3H,0, 
with a:b:c= 0.97636: 1:0.34338, 8 = 89° 18’, has no doubt a closely similar 
structure, with the same tetrahedral framework, but with 2Nat replaced 
by Catt and HO, this replacement lowering the symmetry. The change 














FIGURE 2 


The natrolite framework projected on (001). Natrolite strings, represented by 
squares, are combined by the sharing of unshared tetrahedron corners as shown. ‘The 
number in the center of each square gives the z-codrdinate of the central silicon ion, in 
units c/8, and the numbers at the corners give the z-codrdinates of the shared oxygen 
ions. The heavy lines indicate the longer diagonals of the strings, corresponding to 
the aluminum tetrahedra. 


in composition is of considerable interest—the four cavities occupied by 
2Na* and 2H,0 in natrolite are still occupied, but in this case by 1Ca*+* 
and 3H,0. 

The Scapolites—The minerals of the scapolite group can be considered 
as solid solutions of two end-members, marialite, NasAlsSigOoCl, and 
meionite, CasAlsSigOx4(SOx, COs), in various proportions. They are tet- 
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ragonal, with a:c = 1:0.43925. Data from oscillation and Laue photo- 
graphs of carbonate meionite from Vesuvius and wernerite (of intermediate 
composition) from Bedford County, Ontario, Canada, show that the true 
unit has axes rotated 45° from the crystallographic ones, with 
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FIGURE 3 


The scapolite string. (The distortion 
of the tetrahedra in the upper figure has 
no significance.) 


a = 12.27 A, 
c = 7.66 A. 


The lattice is body-centered, the space 
group being C3,, Cj, or Sj. There 
are two molecules of the above com- 
position in the unit. 

The unit accordingly contains 24 
tetrahedra, of which between 6 and 
12 are aluminum tetrahedra. A rigid 
string of 12 tetrahedra can be built 
as shown in figure 3; if the (Si, Al)-O 
distance is given the average value 
1.65 A, the identity distance along 
the c-axis of the string is found to be 
7.7 A, in good agreement with the 
observed 7.66 A. ‘These strings com- 
bine by sharing still unshared corners 
as shown in figure 4. The structure 
is collapsible. When expanded as far 
as possible it has a = 13.6 A, and, 
moreover, shows vertical symmetry 
planes, which the crystals do not 
possess. Collapse takes place by ro- 
tation of the strings, accompanied by 
loss of the vertical symmetry planes. 
When the angle 0 (figure 4) is about 
25° the edge of the unit becomes 
equal to the observed 12.27 A. 

Assuming all tetrahedra alike, the 
framework is based on the body- 
centered tetragonal lattice, with 
space-group C3,. The space-group 
symmetry may be changed as alumi- 
num tetrahedra are introduced. 

There are two large cavities in the 


unit, occupied by the 2 (Cl, SO,, COs) groups, and there is also room for 


the 8 (Na+, Cat*) ions. 
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TABLE 1 
CALCULATED AND OBSERVED INTENSITIES OF REFLECTION OF SCAPOLITE 
(hkl) ESTIMATED INTENSITY CALCULATED INTENSITY = 42/7120 
0.0.2 1.5 1.57 
0.0.4 50 (50.0) 
0.0.6 10 9.37 
0.0.8 0.5 1.45 
0.0.10 2.0 2.54 
0.0.12 0.0 0.24 
0.0.14 0.2 0.38 








FIGURE 4 


The scapolite framework, projected on (001). Scapolite strings are 
combined by the sharing of unshared tetrahedron corners as shown. 


A first test of the structure is easily made by calculating the intensities 
of reflection from (00 /) with 7 even (the odd orders vanish), using the 
expression 


I~ Abu 
with Ayu = 8Aca + 84s; + 2A4co3 + 16Ao {1 + (— 1)” + cos 2nlz;} 
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+ 16Ag; cos 27lz, in which Al and Si are grouped together. The pre- 
dicted parameter values are z, = 0.175 and z. = 0.205. The agreement 
with the observed relative intensities, obtained by the visual comparison 
of films of varying exposure time, is excellent, and provides strong sub- 
stantiation of the structure (table 1). 

Davynite and Cancrinite—The hexagonal minerals davynite and can- 
crinite are reported by Gossner and Mussgnug‘ to have a = 12.80 A, 





FIGURE 5 


The davynite-cancrinite framework, projected on (001). The circles represent 
oxygen ions at tetrahedron corners; open circles, those with z = 0 and 1/2; solid 
circles, with zg = 1/4; and ringed circles, with z = 3/4. 


c = 5.35 A and a = 12.60 A, c = 5.18 A, respectively, the space group 
being D{,. The davynite unit contains 2 molecules of the composition 
(Na, Ca)sAlSis012(CO3, SOx, Cl). Cancrinite is similar, usually having 
somewhat more silicon than aluminum and containing some water. 

A framework containing 12 tetrahedra in a unit of these dimensions can 
easily be built. The fundamental unit is the ring of six tetrahedra, also 
found in sodalite, beryl, tremolite, and other silicates. If such rings are 
joined by sharing corners as shown in figure 5 there is obtained a hexagonal 
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structure with a = 12.8 A andc = 5.4A. The tetrahedra may be silicon 
and aluminum tetrahedra alternately, each oxygen being common to one 
of each kind. The space-group of this structure is Dj,. There are 
cavities in which the other constituents of the crystal may be placed. 
Indeed, it will be seen that there are tunnels along the c-axis in the form 
of circular cylinders with a diameter of 6 A (after allowing 3.0 A for the 
oxygen ions of the framework). The experimental study of possible 
unusual properties caused by these large tunnels is planned. 

The Cleavage of Tetrahedral Framework Crystals.—The physical properties 
of alkali aluminosilicates are determined in the main by the tetrahedral 
framework of the crystals, for the strength of the electrostatic bond from 
an alkali ion with codérdination number 6 or 8 is '/. or !/s, which is small 
compared with that of an aluminum or silicon ion with coérdination number 
4 (of value */, and 1, respectively). This is strikingly shown by the 
cleavage of these crystals. In table 2.are given the density of AI—O—Si 


TABLE 2 
CLEAVAGE OF TETRAHEDRAL FRAMEWORK CRYSTALS 
NUMBER BONDS PER OBSERVED 
PLANE OF BONDS AREA 100 A? CLEAVAGE 
(110) 4 cv/ 2ab 2.34 Complete 
Natrolite (100), (010) 4 be, ac 3.30 Incomplete? 
(001) 16 ab 4.74 None 
(110) 4 cv/ 2ab 3.01 Good 
Scapolite (100) 4 be 4.26 Less good 
(001) 8 ab 5.31 None 
Davealie: (10.0) 2 ac 2.92 Complete 
Cubic (11.0) 4 acy/ 3 3.36 Less complete or none 
(00.1) 12 a?/3 4.22 Less complete or none 
(110) 4 a?4/2 3.59 More or less distinct 
Sodalite (111) 6 a?+/3 4.39 None 
(100) 4 a? 5.09 None 


bonds (shared tetrahedron corners) which must be broken in order to 
break the framework along the indicated planes, together with the ob- 
served cleavage of the crystals. The correspondence is complete; in 
every case ease of cleavage decreases as bond-density increases. 

A detailed account of the investigation of these minerals will appear in 
the Zeitschrift fir Kristallographie. 


1 Linus Pauling, J. Am. Chem. Soc., 51, 1010 (1929). 

2 Linus Pauling, Zeit. Krist., in press. 

3 Ref. 1, Rule 3. 

4B. Gossner and F. Mussgnug, Zeit. Krist., 73, 52 (1930). 
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THE DEMONSTRATION THAT ALLERGIC INFLAMMATION IS 
NOT NECESSARY FOR THE OPERATION OF ACQUIRED 
IMMUNITY 


By ARNOLD RIcE RICH 


DEPARTMENT OF PATHOLOGY, THE JOHNS HoPpKINS MEDICAL SCHOOL, BALTIMORE, MD. 


Read before the Academy April 28, 1930 


It is well known that, after recovery from many infectious diseases, the 
body which was originally susceptible to infection can be shown to have 
become highly resistant to a second attack of the same infectious agent. 
It is, however, a remarkable fact that this body, resistant to living bacteria, 
is nevertheless extremely sensitive as a rule to noxious effects of products 
of disintegration of the bacteria—so much so that amounts of those prod- 
ucts which would be harmless to the normal body can produce distressing 
symptoms and local damage and death of the tissues, accompanied by 
severe inflammation, in the hypersensitive, or “‘allergic’”’ one. In the case 
of certain infections (tuberculosis, for example) the destruction of tissue 
resulting from this condition of acquired hypersensitiveness is often one 
of the most serious aspects of the disease. It is, nevertheless, commonly 
believed today that the hypersensitive state is necessary for the heightened 
resistance to bacterial invasion which is acquired during infection. The 
idea behind this belief is that since the body is more sensitive to the prod- 
ucts of the infecting bacteria it responds more promptly with inflammation 
on contact with the bacteria, which are thereby prevented from spreading 
and are ultimately destroyed at the site. From this point of view the 
destruction of tissue resulting from allergy is regarded as an unfortunate, 
but necessary local sacrifice for the protection of the body as a whole. 
It is obviously important to know whether the destruction of tissue through 
hypersensitiveness is actually necessary for immunity, for if it is not the 
damage may be diminished or even prevented by methods which are 
effective in reducing or abolishing hypersensitiveness. It is my purpose 
to present in this paper experimental evidence that allergic inflammation 
is not necessary for the operation of acquired immunity. 

Those who believe that hypersensitivity is necessary for immunity base 
their belief upon two arguments. In the first place, it is claimed that the 
two most striking manifestations of acquired immunity—namely, the 
retardation of spread of the bacteria, and the inhibition of their growth— 
can both be explained as effects of the acute inflammation which attends 
the allergic response. This claim is based upon the observation that when 
certain bacteria are injected into tissues previously inflamed by non- 
specific irritants, their spread from the site is retarded and their growth is 
checked. I think that it is important to draw attention briefly to two 
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facts which interfere with our acceptance of the idea that it is the acceler- 
ated acute inflammation of allergy which is responsible for the local fixation 
and the inhibition of growth of bacteria in the immune body. First, 
not all of the types of bacteria which can be restrained by acquired im- 
munity are so affected by a prepared inflammatory exudate. Thus, 
although the activities of the streptococcus, the colon bacillus and bacillus 
pyocyaneus are inhibited when these bacteria are injected into inflamed 
areas, such areas offer no protection whatever against the pneumococcus, ' 
the diphtheria bacillus, the anthrax bacillus,? or the tubercle bacillus.’ 
Second, the injection of bacteria into areas inflamed by the preliminary 
application of an irritant does not actually reproduce the elementary 
conditions of an allergic reaction, although it is widely assumed to do 
so. Inan allergic reaction the bacteria are deposited in uninflamed tissues, 
and a period of time must elapse before the irritant antigen arouses an 
inflammatory response at the site. This situation is comparable to the 
injection of bacteria together with an irritating agent. Using a bacterium 
which is readily inhibited when injected into an area of prepared inflamma- 
tion, Dr. Bull, Miss McKee and I have recently shown that when the time 
relations of allergic inflammation are really reproduced by injecting the 
bacteria at the same time as the inflammatory irritant, in spite of the 
rapid development of inflammation at the site, the growth and spread of 
the bacteria are accelerated rather than retarded, and death of the animal 
occurs earlier than in control animals receiving an injection of bacteria 
alone. As a matter of fact, careful scanning of the protocols of previous 
investigators who have reported that inflammation inhibits the activities 
of bacteria reveals the fact that whenever the conditions of an allergic 
inflammation were reproduced by the injection of the bacteria and the 
irritant at the same time, the promptly developing inflammation afforded 
the animal no protection whatever.”***"’ These considerations prevent 
the acceptance of the claim that it is the prompt inflammation of allergy 
which is primarily responsible for the retardation of spread and the in- 
hibition of growth of bacteria in the immune body. 

The second argument commonly advanced by those who believe that 
allergy is necessary for immunity is that whenever acquired immunity 
is demonstrable, allergy is present also. Since little work has been done 
on the relation between allergy and immunity in infections other than 
tuberculosis, it seemed important to investigate this matter in different 
diseases before accepting the current generalization that allergy and 
immunity are always coexistent. Interestingly enough, in the first disease 
which we chose for investigation, namely, syphilis, we were able to demon- 
strate with ease that acquired immunity can act effectively in the absence 
of allergic inflammation. 

That allergy may develop during the course of syphilitic infection was 
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clearly shown a number of years ago by Noguchi,® and has received ample 
confirmation since. In this disease, as in tuberculosis and in bacterial 
infections in general, the immunity which develops following infection 
has been widely assumed to be the result of allergic hypersensitivity. As 
an example, I need only quote from a recent text book® in which it is 
stated that “immunity (in syphilis) is really a hypersensitiveness.... 
Persons with syphilis are protected from reinfection by the local hyper- 
sensitive reaction which occurs at the site of reinfection.”” In spite of 
the dogmatic character of such statements, as far as I have been able to 
determine there has been no record of an experimental examination of the 
relation of allergy to immunity in syphilis. The study which I am about 
to report was carried out in collaboration with Dr. Chesney and Dr. 
Turner. 

Briefly, we rendered animals immune by the injection of spirochaetes 
derived from human lesions, and at various periods after the immunizing 
inoculation we injected measured doses of virulent spirochaetes into the 
skin of these immunized animals and of normal control animals. The 
results of the injections were studied carefully by gross and microscopic 
examination of the areas at frequent intervals. Under the conditions of 
our experiments, acquired immunity manifested itself not by the assumed 
occurrence of allergic inflammation at the site, but by a remarkable in- 
difference of the tissues to the presence of the injected spirochaetes. 
After a very slight and transient non-specific inflammation which followed 
the injection of the emulsion containing the spirochaetes, and which was 
no differerit in quantity or in quality in the immune animals from what it 
was in the controls, the skin of the animals with acquired immunity re- 
turned to its normal state and remained so, and no late metastatic lesions 
appeared, while in contrast, in the control animals a progressive, ulcerating 
chancre always appeared at the site of inoculation and metastatic lesions 
developed with the usual frequency. The eye and the testis of our im- 
mune animals, sites which are highly reactive in allergic animals, were 
as indifferent to the injection of spirochaetes as was the skin. 

Although, as I have said, allergy undoubtedly may appear under certain 
conditions during syphilitic infection, the present experiments demonstrate 
conclusively, we feel, that in this chronic infection, which resembles 
tuberculosis in so many ways, allergic inflammation is in no sense neces- 
sary for the operation of immunity under experimental conditions in 
animals; and I may add that there is excellent evidence that the same is 
true of human infection. 

What, now, of allergy and immunity in acute infections? For this 
problem the pneumococcus was the organism selected, for it is well known 
that animals actively immunized against the pneumococcus become 
allergic as well as immune, and here, as usual, immunity has been assumed 
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to be dependent upon the hypersensitive inflammatory reaction. I had, 
however, occasionally encountered exceptional animals which had ac- 
quired a high degree of active immunity to the pneumococcus but which 
failed consistently to exhibit allergic reactions on reinfection. These 
accidental observations suggested the possibility of dissociating allergy 
from immunity in this infection by the method of passive immunization. 
Accordingly, Dr. Brown and I introduced the serum from immunized 
animals into the blood stream of normal ones, and then injected measured 
amounts of virulent pneumococci into the skin of these passively im- 
munized animals and of normal, non-immune controls. In no instance 
did an allergic inflammatory reaction occur at the site in the immunized 
animals. All of them survived without developing any appreciable 
lesion at the site of inoculation or elsewhere. In contrast, the non- 
immune controls always developed progressive, local inflammatory lesions, 
and all of them died with septicaemia shortly after the inoculation.!® 

It is clear from these experiments on syphilitic and pneumococcal in- 
fection that it is not an accurate generalization to say that acquired im- 
munity is inseparable from allergic inflammation. On the contrary, the 
experiments demonstrate, in an acute and in a chronic infection, that 
acquired immunity, whether active or passive, can operate effectively in 
the complete absence of allergic inflammation. 

What, then, is the mechanism which retards the spread of bacteria in 
the immune body, and what forces are responsible for their destruction? 
Regarding the first of these problems, I should like to point out that 
bacteria can frequently be found in the regional lymph nodes as early as 
five minutes after having been introduced into the normal body. This 
fact alone should make it clear that the allergic inflammation, which 
requires a much longer time to develop to any appreciable degree, can 
hardly be primarily responsible for the local fixation which occurs in the 
immune body. In a study of pneumococcal immunity I have found that 
local fixation appears to be effected primarily by a prompt and specific 
agglutination of the bacteria which impedes their free movement through 
the tissues; and that this agglutinative fixation (which appears to involve 
a phenomenon of adsorption to the tissues as well) is just as efficacious 
in preventing the spread of the bacteria in passively immunized animals 
in the absence of allergic inflammation as it is in actively immunized, 
allergic animals. It is significant that this process of fixation can occur 
in the tissues of immune animals in which no plasma agglutinins are 
demonstrable by ordinary im vitro tests. Furthermore, while I do not 
mean to suggest that the phagocytes are unimportant in immunity, it is, 
nevertheless, highly interesting that agglutinated pneumococci can be 
seen to undergo destruction and lysis in the tissues extracellularly, without 
the intervention of phagocytosis. 
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In summary, the studies reported in this paper demonstrate that active 
immunity in syphilis does not require allergic inflammation for its success- 
ful operation; that passive immunity to the pneumococcus acts effectively 
in the absence of allergic inflammation; and, finally, that the inhibition 
of spread of pneumococci in the immune body is effected primarily not 
by allergic inflammation, but by a specific agglutinative process which 
acts independently of allergy and is in operation before any inflammation 
appears. 

These facts do not support the current generalization that allergic 
inflammation is necessary for the operation of acquired immunity in 
bacterial infection. On the contrary, they direct attention to the con- 
sideration of abolishing allergy by means of desensitization in diseases, 
such as tuberculosis, in which hypersensitivity is responsible for untoward 
symptoms or tissue destruction. 


1 Clark, Arch. Path., 8, 464 (1929). 

2 Cobbett and Melsome, Cent. allgem. Path. path. Anat., 9, 827 (1898). 
3 Krause and Willis, Am. Rev. Tuberc., 4, 563 (1920). 

4 Opie, Jour. Immunol., 17, 329 (1929). 

§ Rivers and Tillett, J. Exp. Med., 41, 185 (1925). 

® Nakahara, J. Exp. Med., 42, 201 (1925). 

7 Issaeff, Zeits. Hyg., 16, 287 (1894). 

8 Noguchi, J. Exp. Med., 14, 557 (1911). 

9 Low, Anaphylaxis and Sensitization, Edinburgh, 1924. 

1 Rich and Brown, Proc. Soc. Exp. Biol. and Med., 27, 695 (1930). 


ON THE SYSTEMATIC AND ACCIDENTAL ERRORS OF 
MODERN TRIGONOMETRIC PARALLAXES 


By WILLEM J. LUYTEN 
HARVARD COLLEGE OBSERVATORY 
Communicated December 17, 1929 


1. Some years ago (these PROCEEDINGS, 10, 1924, 129-132) E. B. 
Wilson and I discussed two series of parallax determinations, and derived 
a Lexian ratio of 1.24 for the conversion of the published probable errors 
to the actual ones. At that time the amount of parallax material available 
was large enough only in the case of two observatories to be discussed 
in this manner. Since that time, however, large numbers of parallaxes 
have been published by nearly all observatories now engaged in this type 
of work, and it seems possible to apply the same analysis to all photographic 
trigonometric parallaxes now available. 

2. The material used consists of all that could be obtained from pub- 
lished sources by September, 1929, and comprises: 
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Allegheny Observatory (A), Nos. 1-779 with the exception of Nos. 643, 644, 
and 645 which could not be located. 

Royal Observatory, Greenwich (G), all parallaxes from the Greenwich 
parallax volume, and from four lists published subsequently in the Monthly 
Notices, making 411 in all. 

McCormick Observatory (M), Nos. 1-906, with the exception of No. 562, 
which appears to be missing. 

Sproul Observatory (S), Nos. 1-150 in three Sproul Obs. Pub. and an 
additional 22 from Astr. Journal No. 806; 172 in all. 

Mt. Wilson Observatory (W), first to tenth list of parallaxes, numbering 
264. 

Yale Observatory, Johannesburg station (Y), Nos. 1-98. 

Yerkes Observatory (Ye), from all published sources, 292 in all. 

The long list of accurate parallaxes published by the Royal Observatory, 
Cape Town, could unfortunately not be used in the present discussion, 
as they have very few stars in common with other series. 

3. As it is impossible, from relative trigonometric parallaxes alone, 
to derive a systematically correct system of absolute parallaxes, we shall 
confine ourselves here to the calculation of the various systematic differ- 
ences existing between the different series, and to an analysis of the pub- 
lished probable errors and the Lexian ratios to be applied to them. With 
this in view, it is allowable to choose one of the series as standard, arbi- 
trarily. In our case we shall select the McCormick series as standard, 
for the reasons that 

(a) This series is the most numerous. 

(b) It is probably the most homogeneous. 

(c) It has most stars in common with all other series. 

To obtain preliminary information concerning the various types of 
systematic differences that may occur between two parallax series, we 
have first analyzed the difference between the Allegheny and McCormick 
parallaxes, since these, being the most numerous, allow of further sub- 
division, and more detailed analysis. The following sources of systematic 
error have been considered: 

1. Difference depending upon the R. A. of the parallax star. 

2. Difference depending upon the Dec. of the parallax star. 

3. Difference depending upon the apparent magnitude of the parallax 
star. 

4. Difference depending upon the spectral class of the parallax star. 

The results of these comparisons are shown in tables 1 to 4. It might 
be stated at this point that, if a redetermination of the parallax of a star 
should have been undertaken at any observatory for any reason whatever, 
the first value of the parallax has always been used, and the second neg- 
lected. In each table are given, in order, the limiting values of the argu- 
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ment (R. A., Dec., mag., or Spectrum), the algebraic r-~1n values of the 
difference M—A for each interval, expressed in 07001 as unit, and the 
number of stars in each interval. Tables 2 to 4 also contain an estimate 
of the probable error of each mean value, derived simply from the number 
of stars in each interval, and the probable error of one difference M—A 
found from the general distribution of all 329 differences, viz., 0015. 











TABLE 1 
R. A. LIMITS M—A No. R. A. LIMITS M—A NO. 
mhm hmhq*m 
23 30- 0 30 — 7 10 11 30-12 30 + 3 5 
0 30- 1 30 +15 13 12 30-13 30 + 8 5 
1 30- 2 30 +12 10 13 30-14 30 — 4 rf 
2 30- 3 30 — 9 13 14 30-15 30 +12 13 
3 30- 4 30 +8 10 15 30-16 30 + 6 23 
4 30- 5 30 +4 7 16 30-17 30 +1 16 
5 30- 6 30 — 9 13 17 30-18 30 +7 22 
6 30- 7 30 — § 19 18 30-19 30 + 2 26 
7 30- 8 30 — 4 9 19 30-20 30 + 2 32 
8 30- 9 30 +3 15 20 30-21 30 +9 24 
9 30-10 30 —- 9 11 21 30-22 30 +11 15 
10 30-11 30 —12 3 22 30-23 30 — 56 9 
TABLE 2 
DEC. LIMITS M—A No. P.E. 
North of +60° +2.5 20 3.3 
+60° to +40° +1.0 80 ys 
+40° to +20° +2.0 119 1.4 
+20° to 0° +5.1 98 1.5 
South of 0° —0.7 12 4.2 
TABLE 3 TABLE 3A 
MAG. LIMITS M— NO. P. E. MAG. LIMITS M—A No. P. E. 
Br. to 3.0 + 2.2 25 3.0 Br. to 5.0 +1.5 196 ee 
3.04.9 + 1.4 171 1.2 5.0-6.9 +2.6 89 1.6 
5.0-6.9 + 2.0 89 1.6 6.9 and Ft. +6.3 44 2.2 
7.0-8.9 +10 36 2.5 
9.0 and Ft. —11 8 5.3 
TABLE 4 
SPECTRAL CLASS M—A NO. P.E. 
B (+7.0) 2 
A +6.2 46 2.2 
F +3.9 70 1.8 
G +3.5 83 1.6 
K +0.5 109 1.4 
M —5.3 19 3.3 


Viewing the different tables separately, we may make the following 
comments: 
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Table 2.—The ‘ndications are that the systematic difference M—A has 
a tendency to increase from the pole to the equator and to reverse again 
abruptly for southern objects. Considering this reversal, and taking 
into account also the size of the probable errors, it would seem that such 
a general trend in the systematic difference is not yet sufficiently war- 
ranted to give cause for anxiety. It should be added that such a trend 
in the systematic difference is, a priori, not unlikely. 

Tables 3 and 3A.—Owing to the fact that both observatories have used 
rotating sectors to equalize the light of the parallax star to that of the 
comparison stars, a dependency of the systematic difference upon the 
magnitude of the parallax star would appear improbable on the surface. 
Table 3A, which is merely a contracted table 3, indicates, however, that 
the existence of such an error is not imaginary. Its source may possibly 
lie in the fact that for the bright stars the exposure takes place in a different 
way from that for the faint stars, the intervals between each succeeding 
exposure being many times longer than the exposures themselves. It is 
quite possible that a similar correction would be found if the tabulation 
had been made, not on the basis of the original unreduced magnitude 
the parallax star, but on the difference between the reduced magnitude 
and the mean magnitude of the comparison stars. Since these differences 
are small, it requires a much larger amount of data than is now available 
to effect such an analysis with any hope of success. 

Table 4.—Since the Allegheny parallaxes are obtained with a photo- 
graphic refractor, and the McCormick ones through a visual objective 
with a yellow filter, a systematic difference depending upon the color of 
the parallax star is, a priori, not improbable. That such a difference is 
indeed present is plain from table 4, while the size of the probable errors 
attached indicates that this run is not imaginary. Using the effective 
wave-length as a linear color equivalent it might be quite possible to 
establish a workable linear formula to correct for this kind of systematic 
difference. 

Table 1—That a systematic difference depending upon the R. A. of 
the parallax field existed in parallax determinations has long been known, 
and the results of table 1 are here merely given as confirmatory evidence. 
If a sine curve is fitted to the data of table 1, using the numbers in the 
third and sixth columns as weights, the following equation may be ob- 
tained: 


M-—A = +2.0 —2.16 sin a +2.08 cos a 
+0.87 +1.06 +1.35 
giving for the semi-amplitude of the wave 0°0030 + 0”0013 (all errors 


given are probable errors). 
4. Summing up the results of the four tables here analyzed, we might 
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say that a systematic difference between the McCormick and Allegheny 
parallaxes dependent upon spectral class and R. A. of the parallax star is 
definitely indicated, that a systematic difference dependent upon apparent 
magnitude is highly probable, and that one dependent upon declination 
is possible. In view of the fact that investigation of such correlations 
in the differences involving other series is not now feasible because the 
number of stars available for intercomparison is not large enough, we shall 
content ourselves here with simply indicating the existence of these corre- 
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lations in the M—A differences. Correction therefor, before treating 
them together with the other differences, would obviously not be fair to 
the other series, as it would increase the apparent accuracy of the parallaxes 
in the M- and A-series, while leaving the others untouched. It is, more- 
over, impossible, from intercomparison of parallaxes alone, to decide to 
which series the necessary corrections should be made, and it is entirely 
possible that both series should be corrected. 
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From the practical point of view it does not matter much whether the 
M-—A differences are corrected for the systematic errors found above, as 
the resulting decrease in probable error would be small. In first approxi- 
mation, the influence of these systematic errors may be neglected com- 
pared to the much larger accidental errcrs, as is shown by an analysis of 
the frequency curve of all differences. The integrated distribution curve 
of these differences is shown in figure 1, plotted on arithmetic-probability 
paper. Although there are some deviations from the straight line, these 
are so small and erratic, that we may well conclude that the distribution 
of the differences follows the normal law of error. In fact, even at the 
extreme ends of the curve, where differences of —0"047, and +0062 are 
involved, the linearity is well-nigh perfect—sufficient proof, indeed, that 
such large differences should really be expected to occur, and must not 
be rejected, as has sometimes been done by investigators. From the 
straight line we read: lower quartile —0'0135, median +0°0016, upper 
quartile +0°0170. The curve is as symmetrical as could be hoped for, 
and the quartile deviation +0*‘0152 is as nearly as possible equal to the 
probable error found from the squares of the differences, viz., 07015. 

5. In what follows, therefore, the parallaxes have all been used as 
published, without applying corrections of any sort whatsoever. First, 
the five series of Allegheny, McCormick, Sproul, Mount Wilson and 
Yerkes have been discussed; the series of Greenwich, and Yale—Johannes- 
burg requiring special treatment, as will be shown later. The differences 
M —A being already available, the Mount Wilson series was next considered ; 
the reason being that its published probable errors are smaller than those 
of Yerkes and Sproul. The differences M—W, and A—W were obtained 
in the usual way, but wherever determinations from all three observatories 
were available for one and the same star, each difference M— W, and A—W 
was assigned half weight. Subsequently, the Yerkes series was considered, 
and similarly weights of one-half assigned to differences M—Ye, A— Ye, and 
W-— Ye, for stars for which determinations of any three of the observatories 
M, A, W and Ye were available, while a weight of '/; was given when all 
four observatories had determined the parallax of the same star. Finally, 
the S parallaxes were treated in an entirely similar fashion, weights of 
1/, being assigned to each difference involving the S-parallax for stars 
for which all five observatories had made a determination. It has been 
supposed by previous investigators that the S parallaxes show a marked 
dependence upon apparent magnitude; this has been analyzed but no 
appreciable influence found. 

The objection may legitimately be raised that the present procedure 
gives too much weight, relatively, to the series used first, and too little 
to that used last, but it may reasonably be expected that this will not 
influence the results to any marked extent. At least the present method 
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gives the correct total weight to the differences in cases of multiple de- 
terminations, for it would obviously be incorrect to treat all differences 
in such cases as if they were of weight equal to a single difference between 
only two determinations, as has sometimes been done by other investi- 
gators. When five different determinations are available, there are 
evidently only four independent differences, which, from the procedure 
followed here, is exactly equal to the total weight (1 + 2.1/2 + 3.1/3 + 
4.1/4). 

Following the method just outlined, the available knowledge concerning 
differences between parallax series may be condensed into the data col- 
lected in table 5. The first column gives the type of difference, following 
the usual abbreviations for each parallax observatory, the second column 
gives the algebraic mean difference between the parallaxes involved, the 
third column the number of these differences, counted in the way de- 
scribed above, while the last column contains the square of the probable 
error derived from the distribution of these differences themselves (more 
precisely 0.45 of the mean-square of all the differences, uncorrected for 
the fact that the mean value may be different from zero). In the second 
column as well as in the fourth, the numbers are expressed in 0001 as 
unit. 


TABLE 5 
M-A + 2.6 326 222 
M—W — 7.6 62 415 
A —-W — 2.0 29 300 
M — Ye — 2.6 108.3 285 
A — Ye + 0.2 56.9 325 
W — Ye +16.8 21.8 644 
M-S — 9.8 45.5 482 
A -S — 2.5 37.5 216 
w-s — 2.4 9 246 
Ye-S — 4.5 27 464 


Using the numbers in the third column as weights, a least-squares solution 
yields: 


M — A = +1.6 = 1.3 (p. e.) 
M —W = —5.6 + 3.1 
M — Ye = —0.6 = 2.2 
M-S = -—-4.8 = 2.9 


Needless to say, the probable errors of these results have been derived 
by means of the figures in the last column of table 5. Using the mean 
values just derived, we may now correct the numbers in the last column 
of table 5 for the zero point, and thus obtain the real probable errors of 
the distributions of the differences M—A, M—W, etc. According to the 
theory of errors each of these probable errors should then be equal to the 
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square root of the sum of the squares of the probable errors of the con- 
stituents, or, if we denote the probable errors of the various parallax series, 
respectively, by u, a, 8, v, and o, we should get equations of the type: 


pw? + a? = 222-1, etc. 


Again solving these equations by least squares, using the same weights 
as before, we obtain the values for the probable errors a, ...¢ set down 
in table 6. Here the first column contains the designation of each ob- 
servatory, the second the calculated values of a?... and their probable 
errors, the third column the values of a... themselves, with their probable 
errors again. 


TABLE 6 
A 91 + 6 9.5 + 0.3 
M 123 7 Met 0.3 
s 207 20 14.4 0.7 
WwW 262 25 16.2 0.8 
Ye 198 15 14.1 0.5 


Before discussing the significance of these results it would seem well to 
derive, with their aid, the most probable values for the corresponding 
quantities for the parallax series of Greenwich (G) and Yale-Johannes- 
burg (Y). Following the rule given before for weights we find: 

Greenwich parallaxes: 


MEAN VALUE NUMBER SQUARE OF 
DIFF. OF DIFF. OF DIFF. PROBABLE ERROR 
M-—G +3.5 21.5 370 
A-G —5.2 18.9 132 
Ye-—G —5.1 4.8 iets 
W-G —8.4 3.9 


hence, in view of the results from table 5 
M-—G = —1.2; prob. err. of G, 12.6 +1.0 


Yale parallaxes: There are available 29 comparisons with M, 10 with 
A, 2 with W, and 7 with Ye. All the stars in common between Y and W 
or Y and Ye are also in common with M, while of the 10 differences A—Y, 
9 are also in common with M. It would seem simpler, therefore, especially 
in view of the small number of these other three kinds of comparisons, to 
consider only the differences M—Y. For these we find: 


M-—Y = —4.4; square of prob. err., 353 


which would indicate a probable error for the Yale parallaxes of ~/ 230 = 
15.2. (As the number of comparisons is small, it may be pertinent to add 
the corresponding values derived from the still fewer comparisons with 
A and Ye, viz., a mean value of 260, giving a weighted mean of 241 + 45, 
and p.e. = 15.5 + 1.5.)- 
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6. Collecting all results derived thus far, we may combine them to 
form table 7. Here the first column gives the abbreviation of the parallax 
series, the second the systematic correction required to reduce each to the 
McCormick series, and its probable error, the third column the probable 
error actually found for an individual determination in each series, and 
its probable error. 


TABLE 7 
A +0"0016 + 0"0013 0"0095 + 00003 
G —070012 00028 0"0126 0”0010 
M 0” 0000 0%0111 070003 
S —0"0048 070029 00144  0”0007 
WwW —0"0056 070031 0"0162 070008 
Y —0"0044 00035 070155 = 0" 0015 
Ye —0"0006 070022 0”0141 070005 


It is entirely possible that the inclusion of the new Yale parallaxes in 
this investigation is not justifiable, as a systematic difference may well 
exist between parallaxes determined north and south of the equator. It 
would seem too early to comment upon the origin of the systematic differ- 
ences found between the various parallax series, and it may suffice to 
state that we believe those given above to be the most plausible at present. 
Emphasis must be given, however, to their provisional character—which 
to show a mere inspection of the relative size of their probable error will 
be enough. 

In general, we may add that it is gratifying to see that the systematic 
differences between the various series are so small, a fact which speaks 
well of the methods employed in the determination of stellar parallaxes, 
and the extreme care with which measures and reductions are made. It 
is especially reassuring to see that the Greenwich parallaxes appear to be 
systematically only slightly different from the McCormick series. Here- 
tofore the Greenwich parallaxes have usually been omitted by investigators 
discussing parallax series because of the fact that the number of com- 
parisons available was too small to obtain reliable results. Employing 
other means of comparison, large systematic corrections had been de- 
termined such as those of Van Maanen (+0”012 + 070012, Mt. W. Contrib. 
No. 189) and Luyten (+07008 + 07002, Harv. Bull. 799, 1924). Now 
that 49 direct comparisons are available it appears that these indirectly 
derived values are incorrect, and certainly not as reliable as their probable 
errors indicated. 

7. We may also make some remarks concerning the probable errors 
of each parallax series, along the lines of our conclusions in the earlier 
paper. In table 8 a comparison is made between the probable errors 
as published with each parallax (derived from the data given for the 
stars used in this investigation only), the root-mean-squares of which are 
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given in column 2, and the corresponding values actually found here, 
which are taken from table 7 and given in the third column of table 8 
(unit is 0”001). The ratio between them—the Lexian ratio—is contained 
in the fourth column where also its probable error is given. 


TABLE 8 
A 8.0 9.5 1.18 + 0.03 
G 10.0 12.6 1.26 0.10 
M 10.0 EEE LAS 0.03 
SS) 11.8 14.4 1.23 0.06 
WwW 8.2 16.2 1.97 0.11 
= 7.4 15.5 2.00 0.20 
ve 10.8 14.1 1.30 0.04 


Of the two large values, the Mount Wilson one might perhaps have 
been anticipated since these parallaxes were obtained by means of a 
reflector, while of the other, the Yale value, it may again be suspected 
that systematic differences are inherent in parallaxes determined in the 
southern hemisphere as against those observed in the north. When a 
larger number of Yale parallaxes become available, and especially when 
independent data can be derived from the new Cape parallaxes, there 
will be abundant opportunity to investigate this matter further. The 
number of Yale parallaxes used here is admittedly small, but not much 
smaller than the numbers used in the past for comparisons between other 
observatories, before all the modern data were available. 

For the five other series the Lexian ratios, varying from only 1.11 for 
the McCormick parallaxes to 1.30 for the Yerkes parallaxes, are indeed 
smaller than could have been hoped. It may here be recalled that in 
our earlier paper we found, provisionally, Lexian ratio values of 1.24 + 0.06 
for the M—A differences, and 1.4 + 0.08 for the differences involving 
Yerkes parallaxes. With the more extensive material now available 
these have been reduced to 1.18 + 0.03 for Allegheny, 1.11 + 0.03 for 
McCormick, and 1.30 + 0.04 for Yerkes. Considering the fact that the 
published probable errors have been derived from internal discordances 
alone, these Lexian ratios are surprisingly small, thus again providing 
ample proof of the extreme care employed, and the accuracy reached in 
this, one of the most precise kinds of astronomical observation. 
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ON THE DISTANCES AND RADIAL VELOCITIES OF EXTRA- 
GALACTIC NEBULAE, AND THE EXPLANATION 
OF THE LATTER BY THE RELATIVITY 
THEORY OF INERTIA 


By W. DE SITTER 
OBSERVATORY AT LEIDEN 


Communicated May 19, 1930 


1. The large radial velocities of extragalactic nebulae measured by the 
observers at Mount Wilson during the last year have enhanced the im- 
portance of the problem of the determination of the distances of these 
objects. Evidently both the apparent diameter and the integrated ap- 
parent magnitude must be correlated with the distance. In trying to 
derive a mathematical expression for this correlation it is necessary to 
start with a certain number of known distances as a basis. For a few 
of the largest and nearest nebulae the distances have been determined 
directly by Hubble, Lundmark, Shapley and others from cepheids or novae 
in them. For others the brightest stars have been used as a criterion of 
the distance. This criterion, however, is less reliable than the other two, 
because, as Shapley has pointed out, we cannot be sure that what is 
assumed to be a star may not in reality be a cluster of stars. Some of the 
distances determined by this method have, however, been used, as the 
number of distances depending on cepheids or novae alone would be too 
small as a basis for the discussion. Other methods have also been em- 
ployed, based on rotation, or the distance of condensations in spiral arms, 
or other hypothetical considerations; the distances so determined have 
not been used in my discussion. 

2. The first step in the discussion was an investigation of the correla- 
tion between the apparent diameters and integrated magnitudes. The 
discussion is being published in full in B. A. N. 185, to which publication 
the reader must be referred for details. If we assume: Ist, that both 
the true diameters and the true total magnitudes are distributed according 
to a Gaussian (or other symmetrical) frequency function round constant 
average values (constant meaning independent of the distance from us), 
and 2nd, that the measured (or estimated) values of the apparent diameters 
and magnitudes are referred to a consistent scale throughout, then we 
should have, for all objects of one class 


m 
log d + 5 = const., 


where the constant may vary from one class to another. Hubble, in his 
great paper on extragalactic nebulae,’ has assumed this relation. I have 
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FIGURE 1 


Correlation between log d and m/5 for spiral nebulae. The horizontal coérdinate is 
m/5, the vertical coérdinate is log d. Dots are means of from five to fifteen nebulae 
taken in order of increasing m/5, circles are similar means taken in order of decreasing 
log d. The crosses are means of these two kinds of means, and represent the points 
through which the straight line was drawn by Jeast squares. The large circle at the 
top is the Andromeda nebula N. G. C. 224. 
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based my investigation on the data given in that paper which refer to 
291 spirals, 85 elliptical, and 11 irregular nebulae, not counting the “‘un- 
classified”’ or ‘‘peculiar’’ objects. To the irregular nebulae I have added 
the two Magellanic Clouds and N. G. C. 6822, making the total number 
of this class 14. The following formulae were derived. 


m 
Spirals: 3 2.31 — (0.75 + 0.03) (log d — 0.50) 
Elliptical nebulae: e = 2.30 — (0.87 + 0.04) (log b + 0.26) (1) 
1 
Irregular nebulae: = = 1.39 — (0.93 + 0.15) (log d — 1.27) 


e 


For the spirals small corrections have been applied to the diameters and 
magnitudes to reduce them all to the class Sb. For the elliptical nebulae 
the minor diameter b = d (1 — e) is used instead of the major diameter d. 

Although the probable errors in the formulae (1) cannot be assumed to 
be a true measure of the uncertainty, still the deviation from the theoretical 
value 1 is undoubtedly real in the case of the spirals and the elliptical 
nebulae. The data for the spirals are represented in figure 1. It will 
be seen that a line drawn at an angle of 45° would not represent the 
observations satisfactorily. 

The only admissible explanation of the difference of the coefficients 
from unity is by systematic errors in the scale of the diameters, or the 
magnitudes, the small and faint nebulae being estimated too small, or too 
bright, or both. It should be noted that the effect of an absorption of 
light in space would be to make the coefficient larger than unity instead 
of smaller. It is, indeed, hardly to be wondered at that the existing mea- 
sures, or one should rather say estimates, of the diameters and the magni- 
tudes of objects of so widely different aspect as the large and the small 
spirals should not conform to a consistent scale. The problem of estab- 
lishing a standard scale of diameters and of magnitudes for the extra- 
galactic nebulae is a very difficult one, but it is of paramount importance, 
and it is to be hoped that practical astronomers will give it all the atten- 
tion that it deserves. 

3. As has been already stated, the determination of formulas express- 
ing the correlation between distance and apparent diameter, and be- 
tween distance and apparent total magnitude, was based on distances 
determined by Hubble, Lundmark or Shapley from cepheids, novae, or 
the brightest stars. As a convenient unit of distance I use 1074 cm., or 
approximately a million light years, which is denoted by A. The numbers 
of distances used were: 21 spirals (11 Sb and 10 Sc), 6 elliptical nebulae 
(four of which belong to the Virgo cluster, of which the distance is about 
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6A, the two others being N. G. C. 205 and 221 at distance 1A, which 
are, however, not typical for their class) and 6 irregular nebulae. The 
distances of the spirals range from 1A to 6A, those of the irregular nebulae 
from 0.1A to 4.5A. For the spirals and for the irregular nebulae linear 
formulas were derived from these data for the relation between log r 
and log d and between log r and m/5, but the uncertainty of the coefficients 
is naturally very large, owing to the small range of distances. For the 
elliptical nebulae no coefficient could be determined. The two formulas 
for log rz and log 7, were so adjusted as to satisfy exactly the correla- 
tions (1) found between log d and m/5. The finally adopted formulas 
are: 


140 + = 
-140 +2 


Spirals: log ry = 1.28 — 0.75 log d, log ry 


Elliptical: log ry = 0.77 — 0.87 log b, log rm = —1.30 + 3 (2) 


; m 
Irregular: log rz = 1.04 — 0.85 log 6, logr, = —1.30 + 0.91 5 


Thus, for the spirals and the elliptical nebulae, the systematic errors 
of scale, which cause the deviation of the coefficient in (1) from unity, 
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FIGURE 2 

Correlation between distance and magnitude and between distance and diameter. 
The vertical coérdinate is log r. The horizontal codrdinate in the left hand half of the 
figure is m/5, and in the right hand half log d. The broken lines are those derived 
directly from the plotted data. The full lines are the adjusted and finally adopted 
lines, corresponding to the formulae (2). The~diameters of the dots are proportional 
to the square roots of the weights. 








have been thrown entirely on the diameter, the magnitude scale being 
assumed to be correct. This, of course, is rather arbitrary, and probably 
not exact, but no more definite conclusion can be derived from the material 
at present available. The values of log rz and log 7, derived from the 
formulas (2) form a consistent system, of which the scale, however, is 
still very uncertain. The data for the spirals have been represented in 
figure 2. 

The average absolute magnitude according to the formulas (2) is — 15.56 
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for the spirals and — 16.06 for the elliptical nebulae. The former agrees 
well with Hubble’s adopted value. The effect of the spreading in absolute 
magnitude cannot be separated from that of the errors of observation. 
The combination of the two effects corresponds, according to a rough 
determination, to a probable error of +05 for the spirals and +03 for 
the elliptical nebulae. 

4. Comparing the distances thus determined with the observed radial 
velocities, the strong correlation is at once evident. The adopted linear 
formula is: 


log r, = 0.30 = .02 + log (1000 V/c), (3) 


where V is the radial velocity corrected for the rotation of the galaxy, 
and c the velocity of light. For the determination of the coefficient only 
those nebulae were used of which the corrected radial velocity exceeds 
+300 km./sec. Some nebulae, of which the distances or the radial 
velocities were for some reason not reliable, were also not used. Thus 
practically the straight line was drawn through the origin of coérdinates 
and the center of gravity of a group of nebulae at an average distance 
of about 5A. ‘The data are represented in figure 3. The three nebulae 
N. G. C. 4853, 4860, 4865 have been entered at the mean distance corre- 
sponding to the mean of the diameters, the magnitudes being unknown. 
Hubble’s distance for the Coma cluster would bring 4853 and 4860 exactly 
on the curve, while for 4865 the distance derived from the diameter agrees 
with the curve. These three nebulae, of course, have not been used in 
the determination of the coefficient. Also N. G. C. 7619 has not been 
used. 

For the small and faint nebulae there can be no doubt, as has already 
been pointed out by Hubble,’ that the formula (3) gives a better determina- 
tion of the distance than (2). As it depends on (2), it involves the same 
scale, and carries this scale on to large distances, with an error probably 
not exceeding 10%, apart from the errors of observation of the radial 
velocity (which are, however, generally much smaller). 

The meaning of the formula (3) is that the radial velocity is proportional 
to the distance, according to the formula 

SR . 
c 2000" 3) 

The individual velocities are distributed round the value (3’) with a 
spreading corresponding to a probable error of about +140 km./sec. 

5. We call “statical’’ solutions of the field equations of the general 
theory of relativity 


c. are 1/s Suv G +> Zu» + tig = 0, (4) 
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FIGURE 3 


Correlation between distance and radial velocity. The vertical coérdinate is the 
distance, the horizontal codrdinate the radial velocity. Dots are spirals, circles ellipti- 
cal and crosses irregular nebulae. The three dots at distance 40A represent N. G. C. 
4865, 4853, 4860. The two circles near the middle of the figures are N. G. C. 7619 and 
6359, and the cross at distance 16A is N. G. C. 4824. None of these have been used in 
determining the straight line. N. G. C. 3227, 4051 and 6359 are represented in the 
figure, though they were published after the completion of the investigation. 
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those solutions in which the g,, are functions of the three space codrdinates 
Xi, X2, X3 Only, and ‘“‘dynamical”’ solutions those in which they depend 
also on the fourth codrdinate x, = ¢. We restrict ourselves to those 
solutions in which the three-dimensional space has complete spherical 
symmetry, and is filled homogeneously with matter. Under this restric- 
tion the three-dimensional line-element can be written 


R’do? = R*|dx? + sin? x(dy” + sin*ydé*)], (5) 

where 
r gl 
-s (5’) 


is the radius vector expressed in the radius of the universe as a unit, and 
the material energy tensor is 


Tj = — gi3 b, Tu = gu(oo + 39), T = pr, (6) 


where po is the material (invariant) density, which must be a constant 
(i.e., independent of x1, x2, x3), and p is the pressure, which is the sum 
of the radiation pressure of the radiant energy and the dynamical pres- 
sure representing the random motions of the material particles. Although 
in Einstein’s theory the distinction between gravitation and inertia is 
irrelevant, and purely a matter of taste, we can say that, by postulating 
complete homogeneity, we neglect gravitation and consider the field 
of inertia alone. The four-dimensional line-element then is 


ds* = —R*do* + fdt?. (7) 


It is well known that under these restrictions there are only two possible 
statical solutions of the field equations (4) (apart from the limiting case 
x = © corresponding to classical Newtonian mechanics), which I have 
called (A) and (B),? viz: 


1 
(A) R= const = \. f = const. = c’, 


5 (8) 
\" f = c* cos?x. 


In the system (A) the universe is filled with matter, the density depend- 
ing on the radius by the formula 


(B) R = const 


2 
Kpo = R?’ (9) 


where « is the constant of gravitation. The matter has no systematic 
motion:* the path of a test-body (material particle) is a straight line 
(geodesic) described with constant velocity. 
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In the system (B) the density po and the pressure p are both zero:‘ the 
universe is empty. The path of a test-body is in this case a hyperbola, 
of which the asymptotes pass through the origin of coérdinates, and the 
radial component of the velocity is given by 


. # r, v5 

5-pt- S)+9). eee 
7, being the minimum distance, and v, = 17,(d0/cdt)) the transverse 
velocity at this distance.’ Since the time spent by the body in the neigh- 
borhood of its minimum distance is only a short fraction of its whole 
course, we can in all but exceptional cases take as a good approximation 


(B) Pfasaliene (10’) 


Ss 
mI 


We have thus two determinations ef the radius of the three-dimensional 
world, viz.: 


(A) Ri, 


& 
Kpo’ 
(11) 


r 
(B) Rz as V i, c 

The system (A) is only admissible if there are no systematic motions, 
and is thus excluded by the observed systematic radial velocities of the 
extragalactic nebulae. These velocities conform to the law (10’) of the 
system (B), it remaining unexplained, however, why all velocities are 
positive and none negative. This system can, however, only be admitted 
if the density of matter in the universe is so small that p) = 0, or empti- 
ness, can be taken to be a good approximation to the truth. 

Now we know that the part of the universe that we can reach with our 
telescopes is filled with spirals and other extragalactic nebulae. These 
appear to be often condensed in clusters, but as an approximation we may 
take a uniform density of one nebula in a cube of 1A side. Taking the 
average mass of one nebula as unit of mass, the density will be po = 1. 
If for this average mass we take 10" © (estimate by Dr. Oort), and for 
the units of length and of time we take 1A, and the corresponding time 
so as to make c = 1, the constant of gravitation in these units becomes 


x = 3°7.10~7. 


With this value of x, and the value (3’) of V/cr, we find from (11) the 
following values of R, and Rz: 


R, = 2300 A, Rg = 2000 A. (12) 
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The fact that these two values are practically the same, or at least of 
the same order of magnitude, shows that emptiness can not be taken as 
a good approximation in the system (B). The world is, in fact, practically 
full, according to the criterion of system (A). 

We thus come to the conclusion that both the solutions (A) and (B) 
must be rejected, and as these are the only statical solutions of the equa- 
tions (4), the true solution represented in nature must be a dynamical 
solution. 

6. A dynamical solution of the equations (4), with the line-element 
(5), (7) and the material energy tensor (6) is given by Dr. G. Lemaitre 
in a paper® published in 1927, which had unfortunately escaped my notice 
until my attention was called to it by Professor Eddington a few weeks 
ago. In this solution f and R in (7) are assumed to be functions of ¢ 
alone. We can put Vf dt = dr, and use 7 as a new independent variable. 
In other words we can take f = constant = c.? Evidently, if the codrdinate 
x has no systematic motion, which can easily be shown to be the case, 
then the systematic motion of r will be, by (5’), dr/dt = rdR/Rdt, or, 
if by a dot we denote a differential quotient d/cdt, by (11) 


ae ee ‘ 
oR: Re’ (13) 
will be a pure function of ¢, independent of x, x2, x3, of which the value 


at the present moment is, by (12), 0°5.10-* A7. 
The field equations (7) become 


a See 

- ate 

ee See aie 
mt OT, 


, 


where p = po + 3p is the density of “relative mass,’’ while po is the 
density of material, or “invariant,” mass. From these equations we 
easily find 


p+35 ot) =0, (15) 


which is the equation of the conservation of energy. If we introduce the 
volume V = 7°R* of space, this can be written 


d(Vp) + pdV = 0, (15’) 


showing the analogy of the homogeneously distributed matter with a 
gas, from which follows the decrease of the pressure by a diabatic expan- 
sion of the universe. We put further 
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a B 
Koo = Ra» Kp = Rw (16) 


a thus being proportional to the total mass in the universe (a = x«M/z*). 
The equation (15) then becomes 


Ra + 36 = 0. (15”) 


This history of the universe is then described by any two of the equa- 
tions (14) and (15), or (14) and (15’), or (14) and (15”), to which an as- 
sumption regarding a or B must be added. It is most convenient to take 
the second of (14) and (15”). 

7. Lemaitre takes a = constant, from which by (15”) follows B = 
constant. It is, however, certain that material mass is continuously 
converted into energy by the radiation of the stars, and q@ must conse- 
quently have a finite negative value, however small. We can measure 
the rate of conversion of matter into radiant energy against the rate of 
expansion of the universe, and put 

ae 
ae = w R’ 
y being positive. As a convenient hypothesis we can take y to be a 
constant. Then 
a = aR~, 
and from (15”) 
=f ae... Se ae 
B Bo + 3q fu “a 
and consequently 
Bo aoy 1 
= mt 30 — 7) RY 
¥ 





(17) 


The first term of (17) is the kinematic pressure, corresponding to the 
random velocities of the extragalactic nebulae, treated as molecules. 
We have thus for the kinematic pressure 


p _ bo / a _ fo 1 
po sae R Ret a a ROO? 
Now in the kinetic theory of gases p/p is proportional to the square 


of the random velocities. Consequently by the expansion of the universe 
the random velocities decrease proportionally to 1/R'’*~’/*, or, since 








in the actual universe y is a very small quantity, to 1/+/R. 

The second term in (17) is the radiation pressure which is proportional 
to the radiant energy. The total amount of radiant energy in the universe 
is thus proportional to 1/R’, and consequently decreases, i.e., the increase 
of radiant energy by the radiation of the stars is more than balanced by 
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the decrease by the adiabatic expansion of the universe. This can easily 
be verified by making up the account of loss and gain of energy, gain by 
the radiation of matter and loss by the degradation consequent upon the 
increase of wave-length corresponding by Doppler’s principle to the velocity 
of the source. Thus this theory incidentally gives a complete answer to 
the old question what becomes of the energy which is continually being 
poured out into space from the stars. It is more than used up in the work 
done in expanding the universe. It would, however, not be correct to 
say that the universe is expanded by the radiation. It would expand 
just the same if there were no radiation (y = 0). The expansion is due 
to the term with \ in the equation (4), which represents an inherent ex- 
panding tendency of space, counteracting the binding force of gravitation 
represented by x. 

We can make an estimate of the numerical value of y. The total amount 
of energy radiated in the unit of time by a single extragalactic nebula can 
easily be computed from its absolute magnitude. Comparing this with 
the adopted mass we find approximately 


f= —10-1, 


and comparing this with the observed value (13) of R/R, using the value 
(12) of Rg, we have 
y = 2.1077. 

We can thus safely neglect the influence of y, and consequently treat 
a and 8 as constants. 

8. We are now left with the second equation (14), in which po and 
p have the values (16), a and 6 being constants. We put 

=. 
a 
R, being a certain initial value. The equation then becomes 


A fet 2  Z? 
2% 22 + 1 + ale + 22 + b) (18) 





z 


3 2 me ee 
where 

sill inane 

v6 8 —-a—b 

a 6(1+a-— Dd) 

az 3-—a-—) (19) 
B_ (tayo 

R 3-a-)d° 


EE —————————————— 
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The relation between the radius R and the time is given by 


y oa 2dz : 
¥.0- wo f%, (20) 


t) being a constant. Putting 


x=2z-1, 
X? = x? + 4x + B’, B?=3+5), 
the integral becomes é 
(x + 1)dx 
X. V2 + a 


This is an elliptic integral of the third kind. In the case a = 0 it can 
be expressed by logarithms, thus: 


(x + 1)dx 
xX 


: £423 °2 
=Ige +2+X)+ Beery ye 


omitting an additive constant, which can be included in 4. The first 
term becomes positively infinite for x = ©, the second term becomes 
negatively infinite for x = 0, ie, 2 = 1, R = Ro. The radius of the 
world thus increases from the value Rp to infinity, both the initial and 
the final state being reached asymptotically. This is the solution of 
Lemaitre, who, however, only considers the case b = 0, or B = V3. 

An interesting case would be the limiting case B = 2, which gives 
b = 1, a = 0, i.e., a universe devoid of matter but filled with radiation. 
In that case we would have 








(x + 1)dx ———- V(R — Ro)(R + Ro) 
—_—__— = IgV x(x + 2) =1 . : 
wae ev x(x + 2) = Ig z 
This universe would still expand from R = Ry) fort = —~ toR = 
fort = +. 


In the actual universe, however, the value of b is very small. We have 
seen that the random velocities are of the order of 150 km./sec., or 0.5.10~°. 
The ratio p/p = B/aR is proportional to the square of this, from which 
we find, using the actual value z = 2 which will be derived below, 


x = 5.107R, 
0 
or, by (19) 
b(1 + a) = 3.10-° (1 + a — BD), 
or 


b = 3.10 
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The condition a = 0 is the condition that the equation Z? = 0 shall have 
a double root 2 = 1, which must be taken as the lower limit of the in- 
tegral. If a is different from zero, this double root separates into two, 
2 = 1+ ~< —a, which are real if a is negative, and complex if a is posi- 
tive. In the case that a is negative, the lower limit of the integral is the 
largest of these two real roots 2 = 1 + Y-—a,%x% = VY—a; if a is 
positive, the lower limit must be taken % = 0, x) = —1. In both cases 
the integral from the lower limit to an arbitrary value, say x = 1, is 
finite. I find by numerical integration, for the case b = 0, 


1 (x + 1)dx 
01 Vx? — 0.01. X 

1 («x + 1) dx 
~iV x? + 0.01. X 





fora = — 0.01: = ae2 





fora = + 0.01: 3.50. 


By means of the value of \, which will be derived below, the correspond- 
ing values of t — f are found to be 2.3 X 10° years and 4.4 X 10° years, 
respectively, periods which are short in comparison to the lifetime of a 
star. We must thus suppose that in the actual universe the value of a 
is very small. If not exactly zero, it must necessarily be negative, for it 
is inconceivable that the radius of the universe should have started from 
an initial value exactly zero, however long ago that may have been. 

9. From (13) and (18) we have 


a , 
n's* Rs’ (21) 
and from (9), (16) and (19) 

Wite~7} 1 

a Manin - g (22) 


Dividing the square of (21) by (22) we find 
2 


Ra 
Z2=3 757 (lta — d)z, (23) 
Rp 


which is an equation of the fourth degree in z. For b = 0 it becomes of 
the third degree, and if we put also a = 0 it becomes: 


R’, 
23 —_ VA 9 = —. 24 
22— 32+2=3 RB (24) 
The values (12) of R, and Rz give 
Ri 4 
mia 





eee 
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by which the equation (24) becomes 
2 — 3g —2 = 0, (24’) 
from which we find 
z= 2. (25) 


The present value of R is thus twice the initial value Ro. 
Then from (19) and (22) we have 


Meret age 
31 +a— bd) RR 
r= os ) Ra 


3-a-—b 2 
Using the value (12) of Ry, (25) of ganda = 0, db = 0, we find 
d = 1.5.10-* A? 
Ry = 0.8.10 A, 
and consequently the present radius of the universe is 
R = 1.6.10° A = 5.108 parsecs. 


These numerical values are, of course, still very uncertain. 

10. The dynamical solution of the field equations thus is found to 
account for the expansion of the universe, which is observed in the radial 
velocities of the extragalactic nebulae. The universe is homogeneously 
filled with matter and has spherical symmetry throughout its history, its 
radius increasing from an initial value Ry to infinity. The initial state 
corresponds to Einstein’s solution (A). This solution, however, is un- 
stable, and the radius must change. The change is taken to be an ex- 
pansion, in accordance with the observed positive radial velocities, but 
it should be pointed out that the differential equation only gives the 
square of R, and leaves the sign undetermined. ‘This is the fundamental 
problem of why the time has a definite sign, which is inevitably connected 
with the description of nature by differential equations of the second order. 
It should be noted, however, that the dynamical solution requires all 
observed radial velocities to have the same sign, while in the solution 
(B) the sign was indeterminate for each individual body. 

The date of this initial state may, according to the value of the inde- 
terminate constant a, be either a finite or an infinite time ago. The three 
constants 7, a and b have in the theory as here worked out been brought 
into the form of pure numbers, so that their value is independent of the 
units used. It has been found that all three must be very small in the 
actual universe. These constants are of a very different character from 
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the constants A, x, c, which are much more fundamental, and can be 
made equal to unity by an appropriate choice of units. The constant 
X, which is a measure of the inherent expanding force of the universe, is 
still very mysterious, and it is difficult to see what its real meaning is. 
It might even be thought to be one constant too many, unless we may 
hope that it will ultimately be found to be in some way connected with 
Planck’s constant 4. Evidently the dynamical solution of Lemaitre is 
not yet the last word, but it can hardly be doubted that it represents an 
important step:towards the true interpretation of nature. 

LEIDEN, April, 1930 

1 Astroph. J., 64, pp. 321-369; Mt. Wilson Contributions No. 324, 1926. 

2 Proc. Nat. Acad. Sci., 15, p. 172, 1929. 

3 On Einstein’s theory of gravitation and its astronomical consequences, third paper, 
Mon. Nat. R. Astron. Soc., 78, pp. 3-28, 1917. 

4 It should be noted that in both (A) and (B) the pressure is not necessarily zero. In 
(A) there may be radiant energy and random motions of the matter. The complete 
equation (9) is (4p + po) = 2/R*. In (B) the corresponding equation is 4p + po = 0, 
thus if po is positive, the pressure must be negative. See ‘‘third paper,’’ M. N., 78, 
footnote to p. 21. 

5 See “third paper,” M. N., 78, pp. 18-19. 

6 Un Univers homogéne de masse constante et de rayon croissant, rendant compte 
de la vitesse radiale des nébuleuses extra-galactiques, Annales de la Société scientifique 
de Bruxelles, 47, p. 49, April, 1927. A-similar solution had previously been given by 
A. Friedman, Zeitschrift fiir Physik, 10, p. 377, 1922. 


THE CELLULAR DIVISION AND APPROXIMATION OF 
REGULAR SPREADS 


By S. S. Carrns 
DEPARTMENT OF MATHEMATICS, YALE UNIVERSITY 


Communicated June 9, 1930 


The theorems and proofs outlined below are complete and will be elabo- 
rated elsewhere. 

The terms cell simplex and complex are used in the sense of Veblen’s 
definitions. ' 

Let (y) = (1...) ¥n) be rectangular coérdinates in Euclidean n- 
space, m = 2, and let (uw) = (m,..., u;) be codrdinates in an auxiliary 
i-space,1 < mn. A regular i-cell in the space of the y’s will mean an 7-cell 
definable by a correspondence with a simplex in the space of the w’s in 
such a way that 

(1) the y’s as functions of the u’s are continuous with their first partial 
derivatives, and 
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(2) some Jacobian of 7 of the y’s with respect to the u’s is different 
from zero on the cell. 

A closed connected set of points, each of which has for a neighborhood 
on the set a regular 7-cell, will be called a completely regular unbounded 
1-spread. 

THEOREM 1. A region, R, of Euclidean n-space bounded by a finite 
number of distinct completely regular unbounded (n — 1)-spreads is a com- 
plex. 

As a necessary preliminary to dividing R into the n-cells of a complex, 
the points on its boundary are classified according to the direction cosines, 
(a) = (a,...,a@,) of the normal. Class A; (j = 1,...,m) contains 


all points where any of the values | a; | (@@ =1,...,7) is within = of the 


numerical maximum among the a’s, e being a small positive number.’ 
We deal with the points in the order of the classes. 

An n-region of first type extending in the y;-direction will mean an n- 
dimensional rectangle, K},, of the form 


mu < yi < (m; + nj)u 
mu <y¥ < (m; + du (i ¥ j) 


where u is a small positive number and the m’s and n’s are integers, such 
that the bounding (m — 1)-cubes of K}, normal to the y,-axis are free 
from points on the boundary, >>,-1, of R and that K} cuts from )>,~1 
an (n — 1)-cell, single-valued over these bounding (nm — 1)-cubes. The 
bounding 7-rectangles (i = 1,...,m— 1) of Kj}, extending in the y;- 
direction are 1-regions of first type extending in that direction. 

An n-region, Kj, of j'" type (j = 2,.. ., n) extending in the y,-directions 
(«4 = ,..., 4) is defined, inductively, by inequalities of the form 


a <4; a (m; + n;)u (t ™ hi, coe 1;) 
mu < ye < (m, + 1)u (kA t,..., %) 


and has the following properties: 

(1) One face, KZ=}, normal to the yj-axis is an (m — 1)-region of 
(j — 1)th type extending in the y,-directions (¢ = 1, .. ., 4-1), whereas 
the other is free from points on }>,—1. 

(2) Kj is divided by >>, into two n-cells. 

(3) Kj cuts from }>,_1 an (nm — 1)-cell, corresponding by a projection 
to one of the (n — 1) cells into which KZ=} is divided by }>,-1. 

It is also true that any bounding i-rectangle (¢ = 1, ...,m — 1) is an 
i-region of some type from the first through the j“", extending in a subset 
of the y;-directions (4 = %,..., 12;), or else is entirely free from points 


on Binds 
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The y;-direction is possible for an m-region of first type near any 
point of Classes (A;,...,A,). If k directions are anywhere possible, 
then a region of k‘" type extending in those directions is there possible. 
We construct regions of first type in general and use regions of higher types 
to separate those of first type extending in different directions. We thus 
obtain a set of rectangular regions, which contain a neighborhood of the 
boundary of R and each of which cuts a single n-cell from R. The re- 
mainder of R falls into rectangular n-cells entirely in its interior. Bound- 
ing cells are defined so that the totality constitutes an -complex. 

THEOREM 2. A completely regular unbounded n-spread in Euclidean 
(n + r)-space is a complex. 

The points on the spread are classified as in Theorem 1, making use of 
generalized direction cosines of the normal r-plane. We then cut the 
spread into the n-cells of a complex by constructing a set of (m + 1)- 
rectangles of various types, each of which cuts an m-cell from the spread 
and has the further property that any bounding j-rectangle is either free 
from points on the spread or else contains a single bounding cell of the 
n-cell cut out by the (m + r)-region. The defining inequalities of the 
various types of regions are similar to those under Theorem 1. 

A 2-spread will mean a closed connected set of points which fall into the 
two following classes: 

(1) Points each of which has, on the set, a neighborhood which is a 
2-cell. These are called interior points. 

(2) Points, other than interior points, each having, on the set, a 
neighborhood consisting of a 2-cell plus a 1-cell on its boundary. These 
points constitute the boundary of the spread. 

A 2-spread on which the 2-cells of this definition may everywhere be 
chosen so as to be regular, and whose boundary is a 1-circuit made up of 
regular l-cells, will be called a regular 2-piece. A regular 2-spread will 
mean any 2-spread which is the sum of a finite number of regular 2-pieces, 
distinct save for boundary points; and whose boundary, if existent, 
consists of a finite number of 1-circuits, each made up of regular 1-cells. 

THEOREM 3. Any region, R, of Euclidean 3-space, bounded by a finite 
number of distinct regular unbounded 2-spreads, is a complex. 

As in Theorem 1, we divide R into the cells of a complex. 

On any bounding 2-spread, >2, of R, the irregularities are of two classes: 

(Do): Points, each of which is common to three or more regular 2- 
pieces of }°» or else is an irregularity on the boundary of such a 2-piece. 

(D,): Regular 1-cells, each bounded by two of the points Do, but con- 
taining no other points of Do. 

For each point, P, of Do, we construct a 3-cell, cut from R by a cube 
with P near its center. Of the 2-pieces of }>2 and the 1-cells D,, this cube 
meets only those containing P. 





i 
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The irregularities outside these cubes are on a finite set of distinct regular 
1-cells, for any one of which, c,, we construct rectangular regions containing 
a neighborhood thereof. Each region cuts from R a single 3-cell and meets 
only those 2-pieces of }>2 containing ¢. 

After we construct regions near all the boundary irregularities of R, 
it becomes possible to apply the construction methods of Theorem 1 
to the remainder of R and its boundary. 

A regular 2-spread is called completely regular, if every point has a 
neighborhood which is either a regular 2-cell, or a regular 2-cell plus a 
1-cell on its boundary. 

THEOREM 4. For any regular 2-spread Euclidean 3-space and any 
arbitrarily small positive distance, d, there exists an approximating com- 
pletely regular 2-spread in the same space, of such sort that 

(A) The approximating spread coincides with the given spread save at 
points within distance d of the irregularities on the latter; 

(B) A continuous one-to-one correspondence can be set up between the 
two spreads, so that the distance between any pair of corresponding points 
shall be less than d. 

The definition of Do under the preceding theorem applies to interior 
points of any regular 2-spread. We include in this class boundary points 
belonging to more than one 2-piece. The definition of D, applies without 
change to the 2-spread of this theorem. 

The simple case of a 2-spread, }¢s, all of whose irregularities Dy are 
boundary points belonging to just two 2-pieces, is first considered. For 
a 1-cell, c, of D,, a correspondence, of Class C’ save on ¢, is set up between 
a neighborhood of ¢; in 3-space and a neighborhood of some line-segment, 
s1, in such a way that c, corresponds to s,, and neighborhoods of c; on the 
2-pieces of }2 containing it correspond, respectively, to neighborhoods of 
s; on two mutually perpendicular planes. Near s; a regular 2-cell is de- 
fined, whose image under the correspondence replaces a neighborhood 
of c,on }}2. An approximating 2-spread results from the totality of such 
replacements for all the 1-cells Dy. 

The above simple case may be obtained from any regular 2-spread by 
subtracting certain neighborhoods of points belonging to three or more 
2-pieces thereof. An approximation to the simple case thus obtained is 
then extended across its boundary near the subtracted neighborhoods to 
give an approximating 2-spread for the general case. 


1 Oswald Veblen, ‘‘Analysis Situs,” The Cambridge Colloquium, Part 2, page 73. 
2 The idea for such a classification is due to Professor Marston Morse. 
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NORMAL COORDINATES IN THE GEOMETRY OF PATHS 
By Harry Levy . 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF ILLINOIS 
Communicated May 28, 1930 
1. Riemann! first proved that in a metric space characterized by the 
linear element? ¥ 
ds? = g,dx'dx’ (1.1) 
it is always possible to choose a coérdinate system first so that the first s 


derivatives of the g’s are zero at any preassigned point and second so 
that the geodesics through that point have linear equations. Fermi? 
extended Riemann’s results, establishing the existence of a codrdinate 
system in which the g’s are stationary along any preassigned curve, but 
in such a coérdinate system the geodesics will be given by power series ] 
in the arc s in which the coefficients of only the quadratic terms are zero. f 
Veblen‘ and Eisenhart® extended Riemann’s and Fermi’s results, respec- } 
tively, to non-metric characterized by the paths t 
i k 
wha pi, de dx = 0. (1.2) 
ds? ds ds 





Eisenhart’s work has been carried somewhat further by Whitehead and 
Williams. But aside from Veblen no one has found any new codérdinate 
systems in which is retained the vital property of Riemann’s results, 
namely, that geodesics have linear equations. We propose to do this. 
The main purpose of this paper is to prove the following theorem: 

Let V,, be an arbitrary m-dimensional manifold in the space character- 
ized by equations (1.2) and let Ney’ (¢ = 1, 2, ... m) be the components 
ofn —m(o=m+1,m-+ 2, ... nm) arbitrary directions defined at points 
of V,, but none of which lie in V,,; then there exist (infinitely many) co- 
ordinate systems such that the paths through the points of Vm in directions 
linearly dependent on the n — m preassigned directions i.) have linear 
equations. 

If the space (1.2) is Riemannian, the following somewhat special form 
of the theorem seems most useful: 

In any Riemannian space there exist codrdinate systems in which the 
geodesics through the points of a preassigned subspace in the directions orthogo- 
nal to that subspace have linear equations. 

We observe that this theorem may, for m = 0 be regarded as identical 
with Riemann’s. 

2. The proof is direct. Let V,, be given parametrically by the equa- 
tions 
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x = xi(ul, u2, ... u™) $=1,3,...%; (2.1) 


and let \’ be any direction linearly dependent on the ,,)’s? 


¥ = A° Xe)’ (2.2 
The path through a point x in a direction \ is given by* 
Xt = xi + Ws — Yer igd dis? — 1/5, Digs N's3 — ... (2.3) 
where the I'’s are evaluated at the point x. We write 
u’ = A’s c=m+l1,...n, (2.4) 
so that? 
Ns = UW" Ne)® 
and 
ig = x + Uw) = 1/s Tip dcoy Mee) Uw = (2.5) 
If we regard the x’s in these equations as the functions of u', u?, ... u” 


given by (2.1) the X’s are then defined as independent functions of the 
n u’s with single valued inverses in the neighborhood of V,, and hence 
the u’s may be regarded as a new set of coérdinates. 
The path (2.3) determined by the initial conditions 

x = (ui, 2, ... u™) 

rv’ = A’ (uo) Ne) (uo) 
that is, by a point (u}, v3, ... uv’) of V, and a direction there, must sat- 
isfy the conditions 


u* = k* 
uv = k’s (2.6) 
where the k’s are constants. 
We shall speak of coérdinates of this type as normal with respect to Vn. 


3. Let the codrdinates (u!, u?, ... u") be normal with respect to the 
Vn given by 
ue =0 c=m+i1,...n, (3.1) 


and let us introduce new coérdinates a“ defined by the equations 
z= €%(u}, u?, ... u”) (3.2) 


= aru" (3.3) 


n* 
where the a’s are constants and the i’s are functionally independent. 


It follows that the a@’s are likewise normal with respect to V,,. With 
the aid of (3.2) and (2.5) we note that when we make a general trans- 











494 MATHEMATICS: H. LEVY Proc. N. A. S. 


formation of the parameters u!...u” of V,, the normal coérdinates will 
undergo the transformations given by (3.2) and 


“=U o=m-+1,... 2, 


while if we make a general transformation of the space codrdinates or if 
we replace the n — m directions ,,)' associated with the normal coérdinates 
by a new set of )’s the transformation of the normal coérdinates will be 
given by (3.3) and a* = u*, (a = 1, 2, ..., m). 

4. Let Li, be the coefficients of the linear connection determined by 
(1.2) in terms of coérdinates normal with respect to a given V,,. Then 
(2.6) must be a solution of 

dub, 7; du’ du* 
ds? + be ds ds . wad 





whence it follows that 
Li ,k°k" = 0. (4.2) 
Multiplying by s? we find by means (2.6) and the equations immediately 
preceding that the relations 
Liuu =0 i=1,2,...,0 (4.3) 


hold identically throughout the space. Likewise from (4.2), since the 
k’s are arbitrary, it follows that 


pe i=1,2,...,”, (4.4) 
Le, = 0 or =m +1,...9%, 
for all points of V,,. 
If we differentiate (4.3) partially with regard to u” we obtain 
OL, 


T 


uu’ + 2Liu’ = 0, 





whence it follows that the relations 
ais. 


uuu" = 0 4.5 
> (4.5) 





hold throughout the space; again, since the k’s in (2.6) are arbitrary that 
the relations 
p ler _ 9, (4.6) 
ou” 





where P stands for the sum obtained by permuting the indices 7, o, 7 
cyclically, hold for all points of V,,.. We can continue this process and 
obtain a sequency of identities of the types of (4.3) and (4.5) which for 
m = 0 reduce to those obtained by Veblen.‘ 
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5. Let us suppose our space is Riemannian with linear element (1.1), 
that in it we are given a non-minimal subspace V,, defined by (3.1) and 
that the linear element of V,, is given by? 


ds? = Z,edu“du', 


where a bar, as over the g, denotes that the function barred is evaluated 
on V,,, i.e., for uv’ = 0 (o = m+ 1,...n). Weassume, moreover, that 
the codrdinates u' are normal with respect to Vm and that the \’s associated 
with normal coérdinates have been taken as a set of independent mutually 
orthogonal normals. Since the curves of parameters u’ and u", respectively, 
through a point of V,, are orthogonal, it follows that 


age c=m-+l1,...,” 
eae Coe 6.) 


Moreover we may, by a proper choice of u’, make 


Zee = 1 =|,. (5.2) 
If we define the quantities Q(4)ag 
l, og B 
Q ces = — ph pes! 5.3 
(o)as os ae (5.3) 


it follows by a direct computation using (5.1), (5.2) and (4.4) that these 
Q’s are the coefficients of the second fundamental forms of Vi. A similar 
computation shows that the functions p(,¢), defined by 


1 (Gar “Ee 
eas -S par ape : 5.4 
— 2 (S ou" = 





are the coefficients of the linear forms associated by Voss with a subspace.’ 
This interpretation of the 0’s is identical with Bianchi’s form = n — 1,8 
and one could by following his method for that case obtain the Gauss- 
Codazzi equations for a general subspace directly by means of (5.3) 
and (5.4). 

6. Let m = 1 so that V, isacurve C. We assume that space is Rie- 
mannian referred to codérdinates normal with respect to C, and that the 
d’s of §1 are the m — 1 principal normals to C.° The Frenet equations 
for C® 

Jj 
Ne)” J nod ee ee fect Ne-1), 4 fects Ne+1) (6.1) 
ds Po—1 Po 
where the 1/p’s are the curvatures of C, reduce in this case by means of 
(4.4), (5.1) and (5.2) to 
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Ish = et it = — fo, (6.2) 
Po Po -i 
so that 
1 1 
pve t= a 2a) (6.4) 
De 2N Ou Ou" 
or 
1 
— = Ket+1e)1 (6.3) 
Do 


that is, the curvatures of a curve correspond to the y’s of a subspace. 
The analogues of the {’s are, of course, zero. 


1 Gesammelte Werke, 1876, p. 261. 

2 The Latin indices h, i, 7,. .. range through the values 1, 2, . . .m, the Greek a, 8, ¥,... 
through 1, 2, ...m, and the Greek z, o, 7, ... through m + 1,m+2,...n. Repeti- 
tion of an index indicates the sum obtained by allowing that index to take on all values 
of its range. 

3 Rend. Lincet, Rome, 311, 21, 51 (1922). 

4 These PROCEEDINGS, 8, 192-197 (1922). 

5 “Non-Riemannian Geometry,’’ Amer. Math. Soc. Collog. Publ. (1927), p. 64. 

6 Ann. Math., 312, 151 (1930). 

7L. P. Eisenhart, Riemannian Geometry, Princeton, pp. 159-163. 

8 “Tezioni di Geometria Differenziale,’’ Bologna, 2, 450-455 (1924). 

9 Riemannian Geometry, pp. 106-107. 


THE CRYSTAL CLOCK 
By W. A. MARRISON 


BELL TELEPHONE LABORATORIES, NEw YorK CITY 


Read before the Academy, April 29, 1930 


The crystal clock is a relatively new device for keeping accurate time. 
It consists essentially of a generator of constant frequency controlled by a 
resonator made of quartz crystal, with suitable means for producing 
continuous rotation controlled by/it to operate time indicating and related 
mechanisms. 

A crystal clock of this sort has been set up in the Bell Telephone Lab- 
oratories and has been operating over a considerable period. The appara- 
tus was designed especially as a reference standard of frequency for the 
Bell System, but it was recognized that it might also serve as a reference 
standard of time, if developed for that purpose. As a matter of fact, 
since time interval and frequency are so closely related, it would be very 
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desirable to have the same device serve as both a standard of time and a 
standard of frequency. Some types of time measurement can be made 
more readily by means of a time piece of the crystal clock type than by 
any other clock mechanism. This is the chief reason for describing the 
crystal clock at this time. 

The apparatus for a crystal clock consists of a resonator made from 
quartz crystal, with means for controlling its temperature and the sur- 
rounding atmospheric pressure; a vacuum tube oscillating circuit which is 
controlled by the crystal; an electrical circuit, known as a submultiple 
generator, which is used to obtain a low frequency at an exact submultiple 
of the crystal frequency; and a time-indicating mechanism operated from 
this low frequency by means of a small synchronous motor. 

The crystal used is especially designed to have a low temperature 
coefficient of frequency. It is made in the form of a ring with the plane 
of the ring parallel to the optic and electric axes. By using the proper 
proportions, it is possible to make the temperature coefficient of the ring 
as near to zero as desired at a given temperature. This is because the 
temperature coefficient for different modes of vibration is in some cases 
positive and in others negative. By adjusting the shape of the crystal 
the vibration can be apportioned between these modes so that their 
temperature coefficients annul each other. Other shapes® can be used 
to obtain a low temperature coefficient but the ring shape has some addi- 
tional advantages concerned with the mounting. One of the ring-shaped 
crystals adjusted for a frequency of 100,000 cycles is shown in figure 1. 

The crystal is mounted within a temperature controlled chamber and 
kept within about 0.01°C. of a constant temperature. The temperature 
coefficient of the crystal is adjusted to be less than one part in 10° per degree 
Centigrade. The change of frequency due to temperature variations 
alone is, therefore, less than one part in 10°. The crystal with its mounting, 
including the temperature control equipment, is mounted under a her- 
metically sealed bell jar in order to keep the pressure surrounding the 
crystal at a constant value. The effect of pressure on the frequency is 
about one part in 10’ per cm. of mercury. 

The electrical circuit of the oscillator controlled by the crystal, shown 
in figure 2, is somewhat similar to one first described by Professor G. W. 
Pierce.’ The two electrodes of the crystal are connected between the 
grid and filament of the oscillator tube, a grid leak being connected across 
the crystal to provide the proper grid bias. A parallel resonant circuit 
in the plate lead is tuned to a frequency approximately that of the crystal. 
It has been found possible to choose values for the inductance and capacity 
such that a small change in either has a negligible effect upon the frequency. 
Under this condition ordinary changes in filament and plate voltage also 
have very small effects. The oscillating tube is coupled to the output 
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circuits through an intermediate amplifier stage in order to avoid reaction 
on the crystal due to variations in the load circuit. 

A complete crystal controlled oscillator with temperature and pressure 
controlled crystal is shown in figure 3. The dial on the front is used to 
adjust the frequency by means of a small variable condenser connected 
in parallel with the crystal electrodes. Predicted changes in rate as small 
as one part in a hundred million can be made readily by this means. 
The meters indicate the filament and plate currents in the oscillator 
circuit. 

A submultiple generator circuit, controlled by current at the crystal 
frequency, delivers current at a much lower frequency, an exact sub- 
multiple of it. This circuit serves like a reducing gear, in dividing the 
original high frequency by a definite integral amount. In the apparatus 
used at present, the reduction is made in two steps of 10 each, so that one 
cycle in the output corresponds to exactly 100 vibrations of the crystal. 


a 





Ik-4 


+A 
FIGURE 2 


Electrical circuit of crystal oscillator. 


The frequency thus obtained is 1000 cycles per second. This is used to 
operate a clock by means of a small synchronous motor. 

The 1000-cycle motor and clock are shown in figure 4. The rotor of the 
motor is shown in figure 5. The motor unit shown, in addition to the 1000- 
cycle motor, includes a small induction starting motor, two generators 
for producing current at 100 cycles and 10 cycles, and a mercury damped 
flywheel used to reduce hunting. As a clock, all that is required is the 
synchronous motor, some means for starting and a suitably damped 
flywheel. A cam-operated contact device for giving seconds pulses is also 
shown. This is not essential in a clock but is conVenient for making 
accurate time comparisons, as, for example, when checking against time 
signals. 

The frequency reduction in the submultiple generator circuit, and the 
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FIGURE 1 
100,000-cycle low temperature coefficient crystal. 

















FIGURE 5 
_Rotor of 1000-cycle motor. 
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gear reduction in the clock mechanism, are so chosen that, when the crystal 
has its nominal frequency exactly, the clock keeps accurate time. The 
time-keeping properties of the crystal clock, therefore, depend entirely 
upon the performance of the crystal. These, in turn, depend to some 
extent on certain external conditions such as pressure, temperature 
and vibration, which in most cases may be controlled accurately. 

All of the tests made so far on the crystal clock have been over relatively 
short periods of time and it has not been possible, so far, to obtain what 
is perhaps the most interesting data from the standpoint of precision 
clock performance, that is, the undisturbed long time performance. How- 
ever, certain results have been obtained which indicate that there may 
be many uses for a clock of this type and the prospects for further im- 
provement are very good, inasmuch as the total development up to date 
has covered a period of only a few years.”® 

Some performance data indicating the constancy of rate are shown 
in figures 6, 7 and 8. Figure 6 gives the daily gain of the crystal clock 


GAIN IN ONE DAY 
IN SECONDS 
2°29 0 
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FIGURE 6 
Daily rate of crystal clock against Naval Observatory time signals (corrected). 


during December, 1929, and part of January, measured in terms of U. S. 
Naval Observatory time signals. During this period the total change 
in rate was 0.14 second per day, the average variation being very much 
less. ‘Toward the end of the record the observed variations were in the 
order of only one or two-hundredths of a second a day. This record shows 
a gradual increase in the rate of the crystal which has been going on for 
several months. It has not been determined as yet whether this change 
in rate is due to the crystal itself or to some variable in the mounting or 
in the electrical circuit. 

Figure 7 shows the comparison between the crystal used in obtaining 
the data of figure 7 and two others maintained with the same care at 
nearly the same frequency. It can be seen that during this period the 
three crystals remained much more constant relative to each other than 
relative to time signals. 

Figure 8 shows a comparison made between two 100,000-cycle crystals 
to indicate the constancy over short time intervals. This comparison 
was made by means of a ‘special spark chronograph’ which measures 
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accurately the period of each beat between two oscillators and gives a 
final accuracy of about one part in 10". In the figure the greatest de- 
viation from the mean was about one part in a hundred million. 

The crystal clock has a number of inherent advantages over other time- 
keeping devices. Some of these are due to properties of the crystal and 
others are due to properties of the system. 

The properties of quartz crystal which makes it especially adaptable 
for this purpose are, chiefly, the low elastic hysteresis, the low temperature 
coefficient when suitably made, the chemical stability, and the resistance 
to deformation. The stability of shape and composition is perhaps the 
most important factor. Being hard and crystalline, and having been in 
its present form for a very long time, it could be expected not to change a 
great deal, except possibly due to relieving of strains after first being 
cut from the original crystal. A slight aging might be expected also, 
due to the continual vibration, but this can be reduced to a small amount 
by limiting the amplitude of vibration, or might be studied so that the 
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Comparison between the crystal recorded in figure 6 and two similar crystals maintained 
under similar conditions. 


effect could be predicted. Thus, from the standpoint of long continued 
performance, quartz crystal seems to be an ideal material to use as the 
controlling element in the clock. 

In actual tests on the crystal clock, a small aging effect was observed, 
continuing over several months, but decreasing steadily in amount. 
From the available data there is as yet no way to determine with certainty 
whether this aging was due to the crystal itself or due to some uncontrolled 
variable in the system. This will be investigated in further tests using 
fresh crystals. 

On account of the small size of the controlling element and the nature 
of the vibration, it is affected to a relatively small extent by ordinary 
earth vibrations and by slight changes in level. Because of this, it is 
possible to set up a crystal clock in places that would not be at all suitable 
for an accurate pendulum clock. Also, on account of the nature of the 
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controlling element, it is possible to operate a number of them near to- 
gether without any observable mutual interference. ‘The reaction 
between the crystal resonator and its support can be made very small so 
that slight changes in the supporting member itself will have a negligible 
effect upon the frequency of the crystal. 

As far as is known, the frequency of the crystal is not affected by gravita- 
tion or magnetic fields, and it can be shielded perfectly from electro- 
static fields. On account of the independence of the effect of gravity, 
the rate of the crystal would not be expected to vary due to a direct attrac- 
tion of the sun or the moon, or due to a change in gravity caused by earth 
tides. The freedom from effects of magnetic and gravitational fields 
makes it possible to operate a crystal clock in the neighborhood of elec- 
trical machinery or large moving masses. 

With the crystal clock, as with no other high precision clock, it is possible 
to make a continuous adjustment of the mechanism in order that the dial, 
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Comparison of two similar crystals over a relatively short interval. 


indicating the time, shall always be correct within very narrow limits. 
This adjustment can be made either by the use of differential gearing 
in the clock mechanism, or by a continuous phase-shifting device in the 
electrical circuit. With a special form of phase-shifter, the alternating 
current operating the synchronous motor may be advanced or retarded 
gradually by any fraction or by any whole number of cycles. Since one 
whole cycle in this part of the circuit corresponds to 0.001 second, this 
provides a very accurate adjustment. 

The rate of the clock also may be adjusted by predetermined amounts, 
without, in any other way, disturbing its time-keeping qualities. This 
can be done by adjusting a small condenser connected in parallel with the 
crystal electrodes. The effect of varying this condenser is to change very 
slightly the effective stiffness of the crystal. Other variables in the 
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circuit might be used for this adjustment, but the change of capacity is 
the simplest and the most easily controlled. 

Any time-keeping installation, set up for the purpose of keeping accurate 
time, should include at least three clocks so that they can be interchecked 
continuously to expose any erratic behavior that may develop. The 
clocks used in such an installation should be as nearly independent of 
each other, and as independent of common disturbing factors, as possible. 
It should be possible to make continuous records of the intercomparisons 
entirely automatically and with a precision at least one order better than 
that required of the clocks as timekeepers. 

All of these factors can be taken care of quite simply with crystal clocks. 
To compare them against each other with high precision, it is only neces- 
sary to count the beats between the separate oscillators, by means of an 
electrical circuit that has been developed for the purpose, and to record 
the beat numbers at regular intervals. Even the recording of these 
numbers can be done automatically. If one of the crystals is used to 
control an actual clock mechanism, the performance of each of the other 
crystals as clocks can be determined readily, and with extreme accuracy, 
from these data. As a matter of fact, the accuracy of intercomparison 
by this method over a few minutes’ time is as great as the absolute accuracy 
with which the rate of a timekeeper can be determined over a year. 
The crystal previously described makes as many complete vibrations 
in five minutes as a seconds pendulum makes in two years. 

The chief uses for a crystal clock would be those in which it is desired 
to obtain a continuous indication of time, accurate timing signals of 
various sorts, or a continuous and very accurately controlled motion. 
In the event, of course, that the accuracy of the crystal vibration can be 
increased by another order over that obtained at present, there would 
be many obvious uses as an accurate timekeeper, entirely apart from 
the other special properties of the system. 

The crystal clock principle could be used in various problems requiring 
a uniform and very accurately controlled motion. Low frequency currents 
controlled by the crystal can be used to operate synchronous motors at a 
wide variety of speeds, depending on the frequency of the crystal, the ratio 
used in the frequency reducing circuits, and the number of poles on the 
motor. A motor so controlled could be used to control a chronograph, 
high speed oscillograph or even the motion of a telescope. 

A special application of considerable value for the crystal clock would 
be in the transmission of time signals by radio or wire, as is done from a 
number of stations throughout the world. This could be done directly 
from the crystal clock motor by means of a cam-operated contact or by a 
light shutter used in conjunction with a photoelectric cell. If the photo- 
electric cell method were used, signals could be transmitted without the 
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use of physical relay contacts anywhere in the system. By the use of 
amplifiers, the light signal could be modulated directly on the carrier wave. 
Using an auxiliary sectored shutter, or the equivalent, the light going 
to the photoelectric cell could be modulated at exactly 1000 cycles per 
second, controlled by the crystal. With the signals thus produced, direct 
time-signal measurements could be made accurate to a thousandth of a 
second. The beginning of a signal would be indicated by the beginning 
of a pulse, the individual milliseconds would be determined by counting 
the cycles of the 1000 cycles’ modulation from the beginning of a pulse. 
This would involve special indicating or recording equipment but would 
present no real difficulty. It would not in any way interfere with the 
ordinary use of the time signals. Special indications every 10 cycles to 
aid in counting could be provided easily by adjusting the slot widths in 
the light chopper disc. Instead of a complete interruption of tone every 
second, a change of amplitude of say 50 per cent could be used. This 
would indicate the seconds pulses just as effectively as in the present 
method, and, in addition, would allow measurements accurate to a milli- 
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Schematic of proposed time signal transmitter. 


second to be made in terms of the radio signal without interpolating. 
Such accuracy in time measurement would be useful in a study of the 
velocity of wave propagation. One suggested arrangement for such a 
time-signal transmitting device is shown in figure 9. 

It would be easy to insert any identifying signal, such as on the half- 
minute and on the minute. Also it would be possible to divide the minute 
into any number of equal parts, such as 59 to 61, for the purpose of ob- 
taining “rhythmic” signals. 

By the same device that is used to correct the indicated time, it is possi- 
ble also to adjust the time-signal transmitting mechanism so that signals 
can always be sent out as accurately as the time is known. By thus 
eliminating the known clock error and most of the errors of mechanism, 
such as might be introduced by relay contacts, it should be possible to 
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transmit very accurate signals. It might even be possible, except in work 
of the very highest precision, to avoid the necessity for publishing a table 
of time-signal corrections. 

As a timekeeper in astronomical observatories, the crystal clock has 
some interesting possibilities. It would be possible to operate two clock 
mechanisms from the same crystal, one keeping mean solar time and the 
other keeping sidereal time. Both mechanisms are similar, but in the 
simplest mechanical arrangement the gear ratios are different. The 
electrical circuit of one of the frequency reducing units also is changed 
by the addition of a device to make a small correction in the frequency 
controlling the sidereal clock. The ratio of the rates can be adjusted 
in this way as accurately as the ratio of the sidereal to the mean solar day 
is known. According to data given in the American Ephemeris this ratio 
is 1.002,737,803,11 in 1930. The following is a simple design for a com- 
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FIGURE 10 
Schematic of combined mean solar and sidereal clock controlled by a single crystal. 


bined sidereal and mean solar clock adapted to indicate sidereal and mean 
solar time and to produce, if desired, electrical signals for sidereal seconds 
and mean solar seconds. A schematic drawing of such a dual clock arrange- 
ment is shown in figure 10. 

The best fraction with three figure numbers to approximate the desired 
ratio is 366:365. These two numbers factor simply into 61 X 6 and 73 X 5. 
If two synchronous motors have 61 and 73 poles, and are geared down in 
the ratios of 6 to 1 and 5 to 1, respectively, the two slow speed shafts will 
rotate at speeds in the ratio (61 X 6):(73 X 5) = 1.002,739,726,03 when 
the motors are operated from the same frequency. These motors may 
then be geared to clocks through the usual clock train. If the operating 
frequency is such that the clock with the 61-pole motor keeps mean solar 
time, it will operate the other clock at a fair approximation to true sidereal 
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time, the error being about one second a week. The frequency required 
for the solar clock is 366 cycles per second exactly. That required for 
the sidereal clock is 365 X 1.002,737,803,11 = (366 — 0.000,701,865) 
cycles per second. 

The necessary correction in the frequency for controlling the sidereal 
clock is most easily made by means of a continuously operated phase 
shifter connected in the input circuit to the motor, which can be used 
to subtract a small amount from the base frequency. This phase shifter 
may be driven through suitable gearing from either clock motor or from a 
separate motor operating from the same source of current. In the follow- 
ing it will be assumed to be driven from the seconds shaft of the mean 
solar clock motor. The frequency used to operate the mean solar clock 
motor must, in the present case, be 366 cycles per mean solar second. 
This is a little too high for the sidereal clock motor, and 0.000,701,865 
of a cycle per second must be subtracted from it in order to make the 
ratio of the rates correct (according to the astronomical tables) to the 
eleventh significant figure. 

A suitable continuous phase shifter can be built in which one complete 
revolution of the driving shaft adds or subtracts one whole cycle. The 
problem then reduces to that of selecting a gear train that will obtain 
0.000,701,865 revolutions per second from a driving shaft speed of one 
revolution per second. This ratio factors into 1:1000 and 701865:10°. 
The first reduction may be obtained by a double worm gear such as 40 X 
25 or 50 X 20. The second, to a sufficiently good approximation, may be 
obtained by a pair of spur gears having 113 and 161 teeth. The actual 
reduction obtained with these gears, used in combination, is 0.000,701,863 
which gives the desired correction accurate to 0.000,000,002 cycles in 366, 
or better than one part in 10''. With the phase shifter and the gear 
arrangement just described, therefore, the actual ratio of rates obtained 
would differ from the postulated true rate by only one second in thirty 
centuries. 

During thirty centuries the true ratio of the sidereal to the mean solar 
day will probably change by more than that amount. From time to time 
a new gear reducing unit would have to be substituted for the one having 
a ratio of 113 to 161. This would present no difficulty since it is always 
possible to obtain a set of gears that will give any ratio, nearly unity, with 
an accuracy of one part in 10°. 

The current at 366 cycles to operate this dual clock can be obtained 
in a number of ways from a crystal-controlled oscillator as previously 
described. If successive reductions in the submultiple generator were 
taken as 10, 6 and 5 the crystal frequency would have to be 109,800 cycles. 

Such a value would not be so convenient for use as a standard of fre- 
quency as an integral power of 10, such as 10°. If this were an important 
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consideration, a dual clock system as above described could be designed 
to operate from a crystal frequency of 100,000 cycles ot any other power 
of 10. 

In one such possible arrangement both clock moturs have 100 poles 
and rotate at one revolution per second. The mean solar clock motor 
operates from current at exactly 100 cycles obtained from a 100,000- 
cycle crystal by means of a three-stage submultiple generator. The 
sidereal clock motor requires current at 100.273,780,311 cycles to give 
the ratio of rates previously considered. This requires a phase shifter 
rotating about once in four seconds used in this case to increase the fre- 
quency driving the sidereal clock. A good approximation to the desired 
ratio can be obtained with two pairs of gears having the ratios 13:37 
and 60:77 which, used in combination, give a ratio of 0.273,780,274. 
The ratio of rates of the two clocks with this arrangement will differ from 
the assumed true ratio by about one second in a century. 

With the dual clock arrangement, it would be possible to make accurate 
time determinations in terms of either sidereal or mean solar time, always 
using the same time-keeping element. For example, star observations 
could be referred directly to the sidereal clock and any necessary correc- 
tions in the rate could be made in terms of that mechanism by the proper 
adjustment of the crystal oscillator. Such corrections, however, would 
also affect the mean solar clock by just the right amount, so that it would 
be unnecessary to make any computations to determine the corrections 
for the indicated mean solar time. 

It might be possible to make observations on the rate of the crystal 
sidereal clock directly by means of a moving cross-hair, traveling across 
the field of vision in the transit instrument, under the control of a crystal 
controlled motor. The rate of travel would be made the same as that of 
the star image and would be adjusted in space so that, when the clock 
keeps accurate time from day to day, the cross-hair falls exactly on the 
star image. Any error in the clock, therefore, would be indicated by a 
displacement of the star image to one side or the other of the moving 
cross-hair. In order to allow a number of time stars to be used in this 
method, a number of cross-hairs could be used on the same moving element, 
preferably a disc, one corresponding to each favorite star. In using sucha 
method, it would be necessary, of course, to allow for precession and other 
known systematic corrections. Real corrections in the rate of the clock, 
and in the absolute setting of the clock, which should be made from time 
to time, could be made either at the crystal oscillator or in the electrical 
circuit between the crystal and the clocks, so that such corrections affect 
both clocks at the same time and in the proper amounts. 

It would thus be possible to combine, in a single system, mean solar and 
sidereal time-indicating mechan‘s1as, means for rating the clocks in terms 
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of time star observations and means for transmitting time and frequency 
signals with the absolute accuracy of the time determinations. 


The following references cover some of the developments leading up to the crystal 
clock. 


1 Eckhardt, E. A., Karcher, J. C., Keiser, M. Physic. Rev. 17, p. 535, (1921). (De- 
scribes valve maintained mechanical resonator.) 

2Cady, W. G. Proc. I. R. E., April, 1922. (Describes valve maintained quartz 
crystal oscillator.) 

3 Pierce, G. W. Proc. Amer. Acad. Arts and Sct., October, 1923. (Describes im- 
proved crystal oscillator.) 

‘Dye, D. W. Proc. Roy. Soc. London, 103A, p. 241, May, 1923. (Describes tuning. 
fork used as a time standard.) 

5 Horton, J. W., Ricker, N. H., Marrison,W.A. A.J. E. E. Trans., 42, p. 730, 1923. 
(More on fork driven clocks.) 

6 Horton, J. W., Marrison, W. A. Proc. I. R. E., February, 1928. (Describes clock 
operated by a crystal.) 

7 Marrison, W. A. Proc. J. R. E. and B. S. T. J., July, 1929. (Describes crystal 
oscillator with low temperature coefficient crystal.) 

8 Lack, F.R. Proc. I. R. E. and B. S. T. J.,July,1929. (Some propérties of quartz 
resonators, notably temperature coefficient.) 


THE BAND SPECTRUM OF OZONE IN THE VISIBLE AND 
PHOTOGRAPHIC INFRA-RED 


By OLIVER R. WULF 
SMITHSONIAN INSTITUTION AND U. S. BUREAU OF CHEMISTRY AND SOILS 
Read before the Academy April 29, 1930 


In view of the fact that the molecule of ozone has held a prominent 
place in the studies of photochemical and thermal reaction rates, a knowl- 
edge of its spectrum and molecular structure is of interest. It is, further- 
more, one of the simplest of triatomic molecules. The gas absorbs in 
the infra-red, rather generally though weakly across the visible, and, as 
is better known, very strongly in the ultra-violet. In view of this, and 
the fact that it exists in the upper atmosphere, ozone is of considerable 
meteorological interest. In all three regions of the spectrum it reduces 
the intensity of the incoming solar radiation, and it was owing to its 
strong ultra-violet absorption, which constitutes the termination of the 
solar spectrum in the ultra-violet, that the gas was discovered to exist in 
the upper atmosphere. The existing data on the visible absorption of 
ozone does not permit one to gain an adequate idea of the structure of the 
spectrum. In 1880, Chappuis! found eleven bands or regions of absorption 
in the visible and Schoene*? added two more soon after. Later, Laden- 
burg and Lehmann’ studied the visible bands in further detail. Colange* 
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has recently published a curve of the absorption coefficient of ozone across 
the visible and this also shows clearly the positions of absorption found 
by earlier investigators. No attempt appears to have been made, however, 
to look for regularities in the bands, nor to connect them with the rest of 
the known spectrum of ozone. 

In accordance with plans to investigate several cases of especially weak 
molecular absorption, a long absorption path was set up at the Fixed 
Nitrogen Research Laboratory consisting of a Pyrex tube 33 meters long 
and 2'/2 cm. inside diameter. Having such a path available within doors 
this seemed preferable to the use of multiple reflection in shorter paths. 
The ozone absorption was studied in this tube using low concentration 
ozone—oxygen mixtures produced by a silent discharge ozonizer. The type 
of study made it seem inexpedient to work with high concentrations of 
ozone for several reasons. 

A carbon arc burning at 8 to 9 amperes was focussed by a large condens- 
ing lens on the plane window of one end of the tube. Through the plane 
window of the other end of the tube a nearly parallel beam of light emerged. 
A Hilger constant deviation glass spectrograph was used at this end to 
obtain the pictures. The incident light was not parallel but reflection 
from the walls at a high angle conserved much of the light. That which 
entered the walls was for the most part trapped due to total reflection at 
the surfaces. Exposures varied from a few seconds to over an hour. 
As the purpose was simply to obtain good photographs of the absorption 
the absolute amount of ozone in the path was at no time accurately de- 
termined. Wratten M panchromatic and Eastman infra-red-sensitive 
plates were used with the neon and mercury spectra as comparisons. 

It might be clearer to speak of diffuse patches of absorption rather 
than of bands. However, the diffuseness will be discussed below and 
attention is drawn at present to the relative placing of the several bands. 
Plate IA shows a composite made up from different photographs having 
different amounts of ozone in the path and different lengths of exposure. 
Some emission from the arc, in particular three sequences of the cyanogen 
bands, appears on the background but does not lead to any serious confusion 
in perceiving the dark diffuse absorption strips stretching from 6100 A 
to about 4300 A. ‘The plate illustrates fairly clearly the nature of the 
absorption and suggests an interpretation. There are two sets of bands, 
one strong and one weak, and a decided tendency toward regular spacing 
among each set. A possible interpretation would be that the strong bands 
are a progression depending upon an electronic transition from the lowest 
vibration state of the normal molecule, in which transition there is little 
change in the force of binding in the degree of freedom giving these vi- 
brational bands. The regular decrease in intensity in these bands is 
in accord with this. The weaker set would then be looked upon as a pro- 
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gression arising from a low vibrational state of the normal molecule, and 
the displacement to the red of each of these bands with respect to its 
neighboring strong band should indicate the energy of this vibration 
frequency. ‘The separation of the bands in the progressions is of the order 
of 1100 wave numbers while the red displacement of the weak progression 
is of the order of 400 wave numbers.® 

Turning to the infra-red measurements of Ladenburg and Lehmann* 
the strong peak found at 9.7 uw, providing it is a fundamental, accounts 
for a frequency of the order of 1100 wave numbers. In accord with the 
assumption that this is a fundamental frequency, is the band found at 
4.7 » which may be taken as the first overtone of that at 9.7 u, and indeed, 
the peak found at 3.5 u% could be taken as the second overtone, the in- 
tensities being qualitatively in accord. The peak at 6 yu is, according to 
Warburg and Leithauser,® due to nitrogen pentoxide. Regarding the low 
vibrational state of the order of 400 wave numbers, measurements have 
not been made in the vicinity of 254. But the peaks found by Ladenburg 
and Lehmann at 11.3 uw and 7.6 w certainly are not, respectively, a funda- 
mental and its first overtone, yet could be taken as first and second over- 
tone of a fundamental at 23. It seems possible, therefore, that the infra- 
red data indicate the existence of the two frequencies required, but the 
interpretation must evidently be offered only tentatively at this time. 

Owing to the highly diffuse character of the bands it is impossible to 
measure them with much accuracy. They have been measured as well 
as possible from enlargements, the strong bands more accurately than the 
weak bands, and the last bands of each progression least accurately. 
However, with the inability to set upon any significant point in each band, 
the plate with wave-length scale affords about as much information as can 
be given. The values given by the formula vy (cm.~') = 16,625 + n’ 
(1099-20 n’) — mn” (435), for the values »” = 0 and 1, agree fairly well 
with the observed positions, though taking the apparent centers of the 
bands as significant points, the bands do not appear completely in accord 
with a simple formula of this type. The measurements may be taken as 
indicating no inconsistency in the interpretation though with the limited 
accuracy of placing the bands they, of course, cannot be considered as con- 
firmation. 

Extrapolation of the above formula is certainly not justified unless it 
be merely to estimate the order of magnitude of the work of dissociation. 
Linear extrapolation gives about 2 volts as the work of dissociation from 
the excited state, and thus, plus the value of the electronic transition of 
2 volts (16,625 cm.~!), yields about 4 volts for the complete process. 
Taking the 'D level in oxygen atom as 1.4 volts above the normal, the 
dissociation of O; into normal O, and normal O requires 1.3 volts, and the 
dissociation from the electronic level of the visible bands would then 
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appear to be into the '2 state of the oxygen molecule at 1.6 volts and the 
1D state of the atom. 

The diffuseness is, however, of importance in itself. Though the resolu- 
tion is not great the bands not only appear diffuse but the absorption is 
rather general all across the visible. The dispersion is small in the red, 
of course, and relatively large in the blue. At 7600 A the beginning of 
the first atmospheric band of oxygen can be seen, Plate IB, and on certain 
of the original plates with some of its structure resolved. The diffuse 
bands suggest predissociation of the ozone molecule. The process O; —> 
Oz + O requires about 30,000 calories. Thus a quantum of wave-length 
9500 A would have energy sufficient to dissociate the ozone molecule in 
accordance with the above process. All the visible spectrum of ozone 
lies to higher energies, and is overlapped by the continuous of the disso- 
ciation process into normal oxygen molecule and normal oxygen atom. 
This is in accord with the assumption of predissociation. If this assump- 
tion is warranted the initial photochemical act in the photochemical de- 
composition of ozone in visible light is the dissociation of the molecule. 

In the deep red, beyond the beginning of the absorption so far men- 
tioned, further absorption was perceived on certain of the plates. The 
region out to 1 » was, therefore, photographed with neocyanin-stained 
plates. This revealed the series of bands shown in Plate IB. In all, 
ten have been found, forming a series of approximately constant difference 
bands. The formula vy (cm.~') = 10,000 + 566.7 n, where n runs from 
0 to 9, describes them fairly well, though the difficulty in determining the 
positions of the bands, especially those farthest in the infra-red, may not 
permit attaching significance to the formula. They appear diffuse like 
the visible bands, though here the smaller resolution makes this less 
certain. The question of impurities arises, in particular of the oxides of 
nitrogen, since the oxygen was not nitrogen-free. In the presence of 
ozone the only oxide that could be present, except the possibility of nitrous 
oxide, is nitrogen pentoxide. A trap immersed in a solid carbon dioxide 
and alcohol mixture was, therefore, inserted between the ozonizer and tube. 
The bands were still found, however, and as well as could be judged, with 
the same intensity. There are reasons for believing that they are not due 
to N2O, but as this is the next gas to be studied, any uncertainty on this 
point will be cleared later. It seems fair to state that the bands are 
probably due to ozone. 

The above interpretation of the visible bands suggests an influence of 
temperature upon the absorption, thus, a weakening of the weak bands 
relative to the strong ones as the, temperature is lowered,’ and by an 
appreciable amount in passing from room temperature to, for instance, 
—55°C. This suggests a direct method of determining the temperature 
of the atmosphere at an altitude of 50 kilometers. Only very indirect 
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methods at present exist, but the matter has been the subject of much 
calculation and difference of opinion. The ozone fortunately appears 
to be located somewhat as a layer in the vicinity of 50 kilometers. Thus, 
quite aside from the above interpretation of the ozone bands, any tem- 
perature effect on the absorption would afford a possible experimental 
means of determining the temperature of the atmosphere at the height 
of the ozone layer. Indeed it would be important to know of any effect, 
if such existed, since it would influence the determinations of the amount 
of ozone that is in the atmosphere. Since it is also important from the 
standpoint of the structure of the molecule and its photochemistry the 
influence of temperature upon the absorption is now being studied. 


1 Chappuis, Ann. de l’ecole normale sup., [2] 11, 137 (1882). Compt. rend., 91, 985 
(1880); 94, 858 (1882). 

2 Schoene, J. Russ. Phys.-Chem. Soc., 16, Part 9, p. 250(Dec. 20, 1884). J. Chem. 
Soc., 48, 713 (1884); Chem. News, 69, 289 (1894). 

3 Ladenburg and Lehmann, Ann. Phys., [4] 21, 305 (1906). Verh. d. deutsch. Phys. 
Ges., 8, 125 (1906). 

4 Colange, J. phys. et le rad., [6] 8, 254 (1927). 

5 In thus arranging the bands it is taken that there are members of the strong pro- 
gression at 6015 A and 5650 A and members of the weak progression at 6175 A and 
5790 A. It will be noted that the band at 5650 A is weaker than the first member at 
6015 A but stronger than the third member at 5335 A. 

6 Warburg and Leithauser, Ann. Phys., [4] 23, 209 (1907). 

7 Chappuis (reference 1) reported a rather surprising influence of temperature upon 
the absorption of ozone which seems to have been somewhat overlooked, and which 
would appear to be something quite different from the effect suggested here. 


ON THE ESTIMATION OF DISTANCES IN A CURVED 
UNIVERSE WITH A NON-STATIC LINE ELEMENT 


By RICHARD C. TOLMAN 
NorMAN BRIDGE LABORATORY OF Puysics, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated June 14, 1930 


$1. Introduction—In two previous articles,"? I have shown that 
there is a possibility of relating the annihilation of matter in the universe 
to the observed red shift in the light from the extra-galactic nebule, 
by ascribing to the universe a line element of the general form 
ait 2 2 2 
as Se 2 
ds PE sea + dy? + dz?) + di*, (1) 
4R? 


where 7 is an abbreviation for V x? + y? + 2? and Risaconstant. This 
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line element evidently corresponds to a spatially curved universe because 
of the term in r?/R?®, and in addition is non-static because of the presence 
of the time in function g(f). 

The quantities x, y and z, occurring in equation (1) are, of course, 
merely spatial coérdinates which can be thought of as laid off in a 
particular way that was convenient for the derivation of the line ele- 
ment in the form given. And the purpose of the present article is to 
exhibit, more fully than before, the relations between the distance r in 
these coérdinates to a given nebula, and the actual observations of angular 
extension and luminosity that would be obtained by the astronomer for 
the purpose of estimating distances. 

This is a matter of some importance since these relations are obviously 
essential for the theoretical interpretation of the actual astronomical 
observations. Moreover, it will be shown on the basis of the above line 
element that we should expect to find a definite relation between observa- 
tions of angular extension and luminosity, which differs for example from 
that which we should expect for stationary objects in flat space-time, and 
since this difference appears to be just within the existing range of possible 
observational detection, the results of the present paper will afford the 
opportunity for a further check on the correctness of the proposed line 
element. 

§2. Transformation of Codrdinates——As a preliminary to the discus- 
sion, it will first be desirable to give certain forms of the line element that 
result from changes to other systems of codrdinates that sometimes prove 
useful. 

Let us first change to “polar’’ coédrdinates by making the simple sub- 
stitutions 


x = resin 6 cos ¢; y = rsin 6 sin ¢; z= rcos 6. (2) 


We can then evidently rewrite the line element (1) in the form 


eb) 


ds? = — —>——a3 = (dr* + de? + 1° sin? 6dg’) + dt’. (3) 
E + 4R° | 


Let us now make a further change in our method of laying off the spatial 
mesh system by putting 


r i 
oe (4) 
1+ aR 


It can then easily be shown that the line element assumes the form*® 


EE __ — 
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dr? 





ds? = —e + rd + 7? sin? ddg?\ + di. (5) 


r 
— 
This form of the line element will prove to be the simplest to use in the 
actual deduction of expressions for angular extension and luminosity. 
Finally, we shall change to a system of ‘‘rectangular” spatial codrdinates 
by going to a space of higher number of dimensions. To do this we make 
use of the transformation equations 


2, = RV 1-—?r?/R? 2 = r sin 0 cos > (6) 


r cos 6. 


2 = rsin 6sin ¢ Zq 


With these new variables it can be shown that the line element assumes 
the very simple form 


ds? = —ée& (dz? + dzy+ dz? + dz?) + dt. (7) 


In using this form of the line element, however, it should not be for- 
gotten that the actual space of the ‘“‘world’’ occupies merely a three-di- 
mensional ‘‘surface’’ in the fout-dimensional manifold corresponding to 
2... By squaring the transformation equations (6) and adding we 
obtain as the equation of this ‘“‘surface,”’ 


a+ 2+ 23 + a = R’, (8) 


and must remember that the physically significant values of 2... 2% 
will be limited by this equation. We shall later find use for this final 
form of the line element in determining the result of transferring the 
origin of spatial coérdinates from one position to another. 

§3. The Motion of Particles and Light Rays.—As a further preliminary 
to our discussion we need certain information as to the motion of particles 
and light rays that would agree with the above line elements. 

Since the original line element was derived on the assumption that 
particles (nebule), which were stationary in the codérdinates used, would 
remain permanently so, and since the transformation equations (2), (4) 
and (6) are independent of the time, it is evident that the nebule can 
also be taken as stationary in the coérdinates used in the new forms of 
the line element given by equations (3), (5) and (7). 

As to the behavior of particles, which are not stationary, we shall be 
especially interested in the motion of light quanta which are traveling 
radially to or from the origin of the coérdinates used in forms (1), (3) 
and (5) of the line element. In general, of course, the motion of any 
particle will be determined by the equations for a geodesic 

d*x, dx, dx, _ 0. (9) 


gett tuned oe 
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Applying this equation to a particle which is moving at a given instant 
in the x-direction along the x-axis in the original system of codrdinates 
corresponding to (1), we evidently obtain for the acceleration in the 
y-direction 


d*y dx \* : a\? 
aa + {11,2} (+) + {44, 2} (“) = 0, (10) 


and introducing the expressions for the Christoffell symbols given by 
equations (19) in the article already mentioned,‘ this becomes 





d*y 10ufdx\? _ 
ds? 2ody =) = 0, (11) 
where for abbreviation we have set 
est) 
= 9 
od SET, a 
L 4R? 


as in the earlier article. And since for a particle on the x-axis the value 
of Ou/Oy is evidently zero, we see that a particle once moving along the 
x-axis will remain permanently thereon. 

Hence, since the x-axis can evidently be chosen so as to coincide with 
any radial direction of interest, we can conclude as is perhaps obvious 
from the spherical symmetry of the line element, that particles moving 
radially to or from the origin will not be deflected from that path. This 
result, which evidently applies in all three sets of codrdinates (1), (3) and 
(5), and to light quanta as well as to material particles, will be of con- 
siderable use in what follows. 

§4. Transfer of Origin of Codrdinates from Nebule to Observer.—One, 
further preliminary result must still be obtained, which will be needed 
in our discussion of the luminosity of the nebule. Consider a system of 
coérdinates S, corresponding to the line element in the form (5), and 
having a given nebula at the origin of codrdinates and an observer located 
at the coérdinate distance 7 = a, which in accordance with the preceding 
section can be taken as a constant. The needed result now consists in 
showing the possibility of transforming to a new system of codrdinates 
S’, of the same form as before, but having the observer at the origin of 
coérdinates and the nebula located at the codrdinate distance 7’ = a, 
the same numerical value as before. 

To prove this result we shall not attempt to obtain the general trans- 
formation equations connecting the two systems of codérdinates S and 
S’, which turn out to be very complicated, but shall adopt a somewhat 
simpler procedure. In the original coérdinate system S we shall take 


______=—_ 
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the codidinates for the nebula and the observer as shown in the following 
table. 











System S r | 6 ¢ 
Nebula o| - _ (13) 
Observer a | 0 | — 





where for simplicity the observer is given the polar angle 9 = 0, since 
the starting points for measuring the angles @ and ¢ can evidently be chosen 
in any arbitrary way that proves convenient. 

Let us now change to a new system of coérdinates S, which corresponds 
to the line element in the form (7). Making use of the transformation 
equations (6), we then easily obtain the following values for the new co- 
ordinates of the nebula and observer. 


System S, 21 ~ | 2 23 24 
Nebula R | 0 0 0 (14) 
Observer | RV1 — a?/R?|0 | 0 | a 

















Let us now consider a still further change of codrdinates (rotation in 
the 2, 2% plane) to a system of codrdinates S,’, which we define by the 
transformation equations, 





z' = 2,cosa+ x%sina &e! = & (15) 
a’ = —2,sina + 2%cosa 23 = 2 
where we shall put 
. a | a? 
sina = 5 and cos a = ¥! a (16) 


Applying these equations we now obtain the following values for the co- 
ordinates of the nebula and observer. 








System S,’ | 21’ wae ey 
Nebula |RVi1-—a?/R*| 0 | 0 | -a_ (17) 
Observer | R aS ee Se 


It will be noticed, however, from the form of the transformation equa- 
tions (15) that this last transformation has been such as to leave the line 
element still in the form (7) and to preserve the relation (8) connecting 
the physically significant points in the 4-space 2’... 2%’. Hence, we 
may now again use transformation equations of the form (6) to get back 
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to a system of polar coérdinates S’ of the original form (5), in which the 
nebula and observer will be found to have the coérdinates 


System S’ | & 4 ¢’ 











7’ 
Nebula | a | « _ (18) 
Observer | 0 | — - 


Comparing tables (13) and (18), we see that we have now actually 
carried out the desired transformation from a system of polar codrdinates 
with the nebula at the origin to one with the observer at the origin, and 
have proved the equality 7 = 7’ = a of the respective radial coérdinates 
of observer and nebula in the two systems. The result is perhaps an 
obvious consequence of the uniform distribution of matter assumed in 
the derivation of our original line element, but is so important for our 
later considerations as to justify the detailed proof that we have given. 

§5. Estimates of Distance from Observations of Angular Extension.— 
We may now consider the estimation of distances from measurements 
of the mean angular extension of the nebulzinacluster. For the purposes 
of the discussion we shall take the line element in the form given by 
equation (5), and shall take the observer at the origin and the nebula 
to be observed as located at the codrdinate distance 7, so that we may 
regard the light from the edges of the nebula as traveling radially in to 
the observer. Furthermore, we shall take /, as the time at which light 
leaves the nebula and /, as the time of its arrival at the origin when the 
observation is made. 

Let us now take dsp as the diameter of the nebula as measured in proper 
coérdinates by a local observer at the time ¢,, and for convenience let us 
consider this diameter as lying in the direction of 60. From the form of 
the line element we can then evidently write as a relation connecting the 
proper diameter of the nebula 6s) with the codrdinate angle 60, the ex- 
pression 


bso = e*'/? 750, (19) 


where g; denotes the value of g(¢) at the time ¢,. Moreover, since we 
have shown above that light will travel radially toward the origin, it 
is also evident that 66 will be the angular extension ascribed to the nebula 
by the observer at the origin. 

To use the above expression, we may compare this value for angular 
extension with that for a standard comparison nebula of the same in- 
trinsic dimensions, which is located at a standard codrdinate distance 7, 
and emits at time ¢, light which reaches the observer at the same time fp 
as that from the nebula under observation. As an expression for the 
proper diameter of this standard nebula we can evidently write 


i '/i_I!i__ —<— 
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550 = ef" 7,50,, (20) 


where g, denotes the value of g(t) at time ¢,. 
Equating the expressions for the proper diameters of the two nebule, 
given by equations (19) and (20), we can then write 


_ e7 8/2 56, 


's o—e/2 6 (21) 


a 3] 


as an expression connecting the codrdinate distance r of the nebula under 
observation with the value 7, for the standard nebula, and the measured 
angular extensions 6@ and 66, for the two nebule. 

This result, however, can now be put in a more convenient form by 
introducing the expression previously obtained® 





&2— 81 
a (22 





which gives the red shift in the light from a nebula in terms of the values 
of g(t) at the time of departure and arrival of the light. Multiplying 
both the numerator and denominator of the right-hand side of equation 
(21) by e’* and introducing equation (22) we thus obtain 


_ _ (1+ 4A/d) 86, 


ce (23) 
(1 + 4d;/As) 58 


r 
where 6\/X is the fractional red shift in the light from the nebula under 
observation and 6\,/A, is that for the comparison nebula. 

Before proceeding, two criticisms of the above treatment must be 
mentioned. In the first place, the argument has been based on a line 
element which ignores the local gravitational field of the nebula under 
examination, and the light which is finally observed can suffer deflection 
in passing through this local field. In the second place, the nebula selected 
for examination and the comparison nebula emit the light which is finally 
observed at different times and no allowance has been made for possible 
difference in their dimensions due to difference in age. In first approxi- 
mation, however, it is not thought that either of these difficulties is neces- 
sarily very serious, especially since in actual practice nebule which appear 
to be in the same evolutionary stage will be compared, with the result 
that they will tend to have the same proper dimensions and to produce 
the same deflections in light passing through their local fields. 

§6. Estimates of Distance from Observations of Luminosity—We now 
turn to the problem of estimating the distances of nebule from measure- 
ments of their luminosity. For the purposes of the discussion we shall 
again use the line element in the form given by equation (5), but in the 
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first instance shall take the nebula under examination as located at the 
origin of coérdinates and the observer as located at the codrdinate distance 
7, in order that we may regard the light as traveling radially out from 
the nebula. 

As before, let us denote by ¢, the time at which light leaves the nebula 
and by f the time at which it arrives at the observer. In addition let 
us take Zp as the rate at which light quanta of the mean energy Ep are 
leaving the nebula, the measurements being made in proper coérdinates 
which are stationary with respect to the nebula. 

We must now calculate the value L which the observer will obtain for 
the luminosity of the nebula, which we define as the energy received from 
the nebula per unit time and per unit cross-section as measured in the 
local coérdinates of the observer. ‘To do this we note, from the form 
of the line element (5), that the proper area of the spherical surface at 
7 through which and perpendicular to which the quanta are passing at 
time ¢ is evidently 47 7? e*; that the average energy of the quanta as 
measured at this surface will evidently be E)/(1 + 6A/d); and finally 
in accordance with our previous work® that a short time interval dt, during 
which quanta leave the nebula will be connected with the time interval 
dt, in which they reach the surface at 7 by the relation 





d, =* 


—mse* =] Od/X. 24 
anne + 8n/ (24) 


With these considerations in mind we can then evidently write for the 
measured luminosity of the nebula 


L a Zo Eo dt, a ZoEo 1 
~ dpe 14+ 5d/X df, 4ar?e® (1 + 6X/d) 


And for the luminosity of a similar standard comparison nebula we may 
write 








(25) 





ZoEo 1 
L, = 26 
© Aw te (1 + 5d;/d;)? (26) 
or combining the two equations we have 
1+ aya) i. 
7 =7, (1 + d;/As) ni (27) 


(1+ 6x/ra) VL- 


In accordance with the method of derivation the quantity 7 occurring 
in this equation locates the observer in a system of codérdinates which was 
taken for convenience with the nebula at its origin. Nevertheless, since 


se 








Vo. 16, 1930 PHYSICS: R. C. TOLMAN 519 


we have already shown in §4 that the radial coérdinate for the nebula 
in a system of codrdinates with the observer at the origin will have this 
same value, we can take equation (27), as locating the nebula in a system 
of codrdinates with the observer at the origin. 

Concerning the validity of equation (27), similar remarks can be made 
with regard to the neglect of local gravitational fields and the possible 
effects of a difference in age of the nebulae as were made concerning equa- 
tion (23). In addition it should be pointed out that the luminosities in 
equation (27) are obviously bolometric luminosities and corrections to 
pass from visual or photographic to bolometric magnitudes must be made 
if possible, before applying the equation to the usual astronomical ob- 
servations. 

§7. Conclusion—The above equations connecting the codrdinate 
distance 7 to the nebula with angular extension 6@ and luminosity Z were 
derived using the line element in the form (5), but with the help of the 
transformation equation (4) they can, of course, also be expressed in terms 
of the codrdinate distance 7 occurring in the line element in the forms 
(1) and (3). For convenience we collect the results below 


__ (+ an/r) _ (1+ 4r/r) IG 
. "Gt Gia e.. Gee Veo 
(1 + r?/4R?) (1+ 5d/d) 460, 


T's C+ /4R®) (1+ 8, /2s) 80 


(1 + 7°/4R*) (1+ 8rs/ds) o|Ls 
sagas 2 2 yes (29) 
(1 + r°/4R3) (1+ 6r/r) VL 
These are the equations to be used when we desire to obtain the coérdinate 
distances of the nebule from the results of actual measurements of angular 
extension, luminosity and red shift. 
Finally, it should be specially pointed out that the equality of the two 
different expressions for coérdinate distance given by either (28) or (29) 
evidently leads to a relation between angular extension and luminosity 


(28) 








6 
= const. X (1 + 65A/d)? (30) 


which could be made the subject of experimental test. Furthermore, 
since the measurements have already been pushed to a distance where 
the value of the red shift 5\/ has risen to 0.04, it is hoped that this method 
of testing the theory may actually be possible in the not too distant future. 
The check on the theory thus provided would perhaps be a minor one, 
and the interpretation of the observations would be subject to the diffi- 
culties referred to at the end of §5 and of §6, but the undertaking would 
certainly be an interesting one since the relation, although agreeing with 


' 
“om 
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that which would be expected for objects in flat space-time with an actual 
motion of recession, differs from that which would be expected on the basis 
of some explanations of the red shift which might be proposed. 

1 Tolman, Proc. Nat. Acad. Sci., 16, 320 (1930). 

2 Tolman, Jbid., 16, 409 (1930). 

8 This is the form more similar to the usual Einstein one mentioned in my previous 
article (Ref. 1, p. 334) and the form of Robertson (see Ref. 2, § 2). 

4 See Ref. 1, p. 325. 

5 See Ref. 2, equation (21). 


EXPERIMENTS IN BIRD MIGRATION. II. REVERSED 
MIGRATION 


By WILLIAM ROWAN 
DEPARTMENT OF ZOOLOGY, UNIVERSITY OF ALBERTA 
Communicated May 23, 1930 


The experiments in bird migration, which culminated during the winter 
1929-1930 in an attempt to induce the common crow, Corvus brachy- 
rhynchos brachyrhynchos Brehm, to migrate north instead of south in the 
fall, began in the autumn of 1924. The bird then used was the Junco, 
Junco hyemalis connectens Coues, a common breeder in central Alberta, 
available in large numbers, hardy and a fairly extensive migrant. The 
“working hypothesis’’ on which the experiments were based may be 
briefly outlined as follows. It was assumed that the migratory habit 
was, in the majority of migrants of the northern hemisphere, inherent 
and that food, cold and all the other factors, all of them highly variable, 
commonly held to account for the southward passage could not account 
for it except possibly in the case of certain races of a few species. Most 
species migrate with great regularity whatever the weather may do. In 
the case of the earlier migrants, leaving the north in July and early August, 
stress of circumstances can certainly not be cited as the stimulus to the 
southern urge. The accuracy from year to year of the dates of departure 
suggested that some environmental factor, which must be stable, was never- 
theless involved as a timing agent. It was argued that this might well 
be found in the autumnal day-lengths. This, however, would not be 
sufficient to induce the passage if the tendency to undertake it were not 
inherent. An internal stimulus, affecting the physiological condition 
of the bird and arousing the impulse to depart also had to be stipulated. 
Internal secretions of the reproductive organs suggested themselves for 
the reason that the gonads of birds show a remarkable cycle coincident 
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with migration. They are minute in the winter and relatively enormous 
in the summer, the bird migrating northward as the organs enlarge and 
southward as they decrease. Most immature and non-breeding birds, 
whose gonads do not reach the full condition, fail also to migrate or achieve 
only a partial migration. The ‘working hypothesis,’’ then, assumed 
that day-length controls the condition of the reproductive organs and 
that these at certain stages produce hormones which in turn afford the 
physiological stimulus that inaugurates the north and south passages. 
Physiologists and popular authors have assumed that the reproductive 
cycle of birds is dependent on temperatures. But this viewpoint, for a 
variety of reasons that need not be detailed, appeared to the writer to be 
quite untenable. When the experiments were undertaken heated aviaries 
were, therefore, not bothered with and in order to put them to the severest 
possible test, the birds were housed out of doors where they could be 
expected periodically to be subjected to temperatures running down to 
50°F. below zero. The experimerital aviary was lighted by means of 
electric light in such a way that, commencing early in October, the birds 
had days increasing in length at the rate of 5 minutes per day while the 
controls, without artificial light, were getting the shortening days normal 
to the fall. Samples were taken at intervals. The gonads of the controls, 
after reaching the absolute minimum in early November, remained there 
till the end of February. Those of the experimentals soon began to 
show an increase in size and by New Year had attained the full spring 
condition in spite of the most drastic temperatures, —52°F. being the 
lowest recorded during the running of the experiments. 

The next step was to liberate controls, with their gonads in a state of 
rest and experimentals, with their organs physiologically active and 
increasing in size and to record their behavior. According to the ordinary 
canons of migration, both lots should have gone south on liberation since 
they were released under winter conditions several hundred miles north 
of their wintering grounds. On the other hand, according to the hy- 
pothesis, the controls should have made no attempt to migrate at all, 
while the experimentals should have gone north, their organs being in an 
active state of recrudescence. Controls liberated from November to 
February proved entirely sedentary. Except for three or four birds almost 
certainly secured by cats, all were retaken in traps at the aviaries. Of 
the experimentals, anything up to 60% left the premises for good on 
release. For full details and discussion reference may be made to Rowan, 
‘Experiments in Bird Migration,” Proc. Boston Soc. Nat. Hist., 39, No. 5. 

For a variety of reasons it was particularly desirable to ascertain, if 
possible, where these experimental birds went on departure. The Junco 
was unsuitable for the purpose on account of its small size, the fact that 
it is protected by law and because few people know it when they see it. 
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It would have been useless, and illegal, to invite people to collect the birds 
after release and return them tous. A large migratory species, well known 
to everyone, unprotected by law and able to find suitable food on libera- 
tion, was obviously needed. The crow seemed designed for the purpose, 
the only objection being the supposed impossibility of securing it alive 
during the summer months. . All efforts in Canada and the States known 
to us have been failures and we were universally recommended not to 
waste time attempting it. (Jack Miner traps crows only during the 
winter months.) However, on August 14, 1929, trapping was started. 
By mid-September, when almost the last crows had gone south, we had 
140 birds, gotten by a considerable variety of methods. This number 
would certainly have been trebled had our main trap not been interfered 
with at the crucial moment when crows were going into it, scores at a 
time, to the lure of rotten eggs. 

Our original plans were as follows. In view of the geographical position 
of Edmonton near the northern limit of cultivated Alberta, with large 
tracts of muskeg and uninhabited forest to the north and northwest, 
it was intended to run the experiment at least 100 miles to the south or 
southeast. It was further planned to castrate a percentage of the ex- 
perimentals and to use injections of testicular emulsion both on birds so 
treated and on some of the controls. All this had to be abandoned on 
account of the small numbers available. Some died on first being caged; 
others were taken as samples from time to time, leaving us finally with 
but 83 crows at the time of liberation on November 9. The interference 
of routine duties would have made supervision of the lighting and general 
care of the birds impossible at some distant point and in the absence of 
a man with previous experience who might have taken control, Edmonton 
was reluctantly decided upon after all. 

The birds were so divided that on Nov. 9 there were 69 experimentals 
and 14 controls. Both groups had a large flying-cage attached to a shed 
for shelter, the former, since they were by far the more numerous, having 
somewhat roomier accommodation, their flying-cage being 100 ft. in length, 
6 in height and 15 in width. Food for both lots was identical and con- 
sisted in the main of rotten eggs, dogs and cats from the city pound, an 
occasional horse, green stuff of various sorts, grain and chicken mash, 
stale bread and buttermilk, etc. This seems to have constituted a well- 
balanced diet. The birds could not have been in better shape at the time 
of liberation. 

The experimentals received lighting on the same basis as their prede- 
cessors, the Juncos, the cage and shed being illuminated by means of 
twenty-five 500-watt bulbs, a total of 12,500 Watt. Lighting was begun at 
sunset (6.15) on September 28 and continued till 7.15. Thereafter the 
lights were turned on daily at sunset and continued for 7!/2 minutes longer 
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than on the previous day, the effective daily increases thus being approxi- 
mately 5 minutes (allowing for the shortening of the fall mornings.) On 
Oct. 20 they were burning till 10 P.M. Thereafter the nightly increase 
was reduced to 5 minutes till Nov. 1, when 11 o’clock had been reached 
beyond which the time was never extended. On the whole the birds 
stayed awake very well although there appeared to be some individual 
variation. Many of them fed by the artificial light as during the day, 
though the frequent arrival of visitors tended to confine them to the shed 
in which no food was available at any time. 

At the close of October the tester showed the first indications of de- 
veloping interstitial tissue and it was then decided to take advantage of 
the large turn-out of gunners likely on Thanksgiving Day, Nov. 11, and 
the prospects of good roads (there had been no snow up to date) and to 
liberate the birds in time to have a first collection made on that day. 
The 9th wag a Saturday. Since shooting is entirely prohibited in the 
Province on Sundays, the birds “were surreptitiously released on the 
morning of the 9th and the fact announced in the evening newspapers 
and over the radios, with an appeal for anyone and everyone to hunt 
crows on Thanksgiving Day (and thereafter) and to mail the birds in to 
the Department of Zodlogy at the University. All birds were banded. 
The crows thus had some 48 hours in which they could disperse and 
start on their respective roads before being hunted. The date was earlier 
than could have been wished, particularly with respect to the controls, 
but it seemed a chance not to be missed. The outcome, however, has 
shown that Thanksgiving Day brought but slight benefit, in spite of a 
truly remarkable response. Hundreds of people were on the lookout for 
crows that day, but the birds seemed to be just as wild and wary as when 
first trapped and impossible to approach. At least 90% of those ulti- 
mately shot were secured by residents in various places onto whose premises 
the birds had strayed. 

As far as the controls were concerned, it was doubtful if the gonads had 
actually reached the minimum by Nov. 9. It seemed essential, however, 
that they should be liberated under exactly the same conditions as the 
experimentals and both groups were therefore liberated simultaneously. 
The flocking habits of crows during the fall and winter was a serious ob- 
jection to this double release, as the behavior of the one group might well 
be expected to affect that of the other, but it seemed to constitute a more 
severe test than separate liberation could have done. If the birds sorted 
themselves out into groups, such segregation would prove the more con- 
vincing. 

In the following summary of the results, birds shot within a ten-mile 
radius of Edmonton are included as local kills. It is assumed that the 
direction taken beyond this limit is significant, but how far the assumption 
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is justified is debatable. The distance should possibly be considerably 
increased and the time factor taken into consideration. 

Of experimentals, 69 were liberated on Nov.9. At9A.m. the barometer 
(rising) stood at 27.30 (which is low for Edmonton) with an official forecast 
of “‘local showers or snow flurries.” It was cold but not freezing, cloudy 
with a light wind. Of these 69, 28 stayed and were retaken on the 11th 
and 12th. The remaining 41 left us for good. A constant watch for 
possible returns was kept throughout the winter. The ones that stayed 
were given lighting as before for an additional two weeks and then turned 
out again under less favorable weather conditions. Thirteen departed 
while the remaining 15 were retaken within two days. Thus, of the 
original 69 birds, 54 (78%) ultimately left us. Of these 54, 28 have been 
reported killed, either directly to us or through the Biological Survey at 
Washington, whose address was stamped on the bands. The remaining 
26 are unaccounted for at time of writing, May 14, 1930. Of the 28 
reported killed, 12 are local (within 10 miles of Edmonton), 8 south 
and southeast (the farthest in South Dakota, U. S. A.) and 8 north 
and northwest (the farthest—-2 birds—at Whitecourt, Alberta, 100 miles 
“as the crow flies’ north of west of Edmonton.) Of our 54 experi- 
mental departures, 8 (15%) have thus been reported from the south and 
southeast. 

Of the 14 controls liberated on Nov. 9, 6 remained and 8 departed. The 
ones that remained were liberated again on the 24th when none left us. 
Of the 8 departures, 6 have been reported killed, 2 locally and 4 to the 
southeast (the farthest about 200 miles). Two are unaccounted for. 
Of this group, 50% of the departures have thus been reported from the 
southeast and none from any point north or northwest. 

Brief comment may be made on the above facts. The total number of 
crows involved was unfortunately small, an unavoidable circumstance. 
In addition to this two other factors militated against decisive results— 
the geographical position of Edmonton and the early date of release. As 
far as the controls are concerned, there is every reason to believe that had 
they been kept another two weeks, all would have stayed. The small 
number of birds kept as controls was due to this expectation. It was also 
eminently desirable to have the largest possible number of birds in the 
experimental group if an adequate idea of the directions taken by them 
was to be obtained. 

If the experimentals had actually gone northwest they had before them 
sparsely settled country from which but few returns could at best be 
looked for. The securing of two birds at Whitecourt, in bush and muskeg 
country, is really remarkable. The first of these was taken on Nov. 16, 
just a week after liberation: the second on February 28, 1930, nearly 
four months after release. The latter is a particularly interesting record 
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since it demonstrated the ability of the bird to fend for itself and to survive 
our drastic winter in the wilds. 

While sight reports sent in of crows have been numerous those from the 
south have, with one exception, referred to single crows or rarely two 
together. From the northwest we have, on the other hand, had reports 
in the main of groups of crows (up to 15 together), reports from the Lesser 
Slave Lake country, as far as Peace River Town, 300 miles northwest of 
Edmonton, having been particularly persistent till January. That there 
were crows up there through most of the winter is beyond doubt. The 
birds were first noted in the locality on the 18th of November, 9 days after 
our releases and nearly two months since the last crow had previously 
been observed here. Several efforts were made to secure samples but 
all failed, chiefly on account of the wildness of the birds and the impedi- 
ment of deep snow. In view of the fact that 50% of the controls were 
recovered in the well-settled country south of Edmonton and only 15% 
of the experimentals, it is difficult to believe that anything but a minority 
of the latter group could have gone south. One is left to infer that the 
Lesser Slave Lake crows represent, at least some of, the missing 48% of 
our experimental birds. It is most regrettable that this particular group 
of chief interest, had they fulfilled expectations, should have been destined 
to head into the wilderness where returns would inevitably be few. But the 
paucity of returns from the south (15%) where an excellent watch was 
kept, coupled with the frequent sight records from the Lesser Slave Lake 
territory, distinctly suggests that they actually did what was anticipated. 

In view of the doubtful situation that now exists and the suggestive 
circumstantial evidence, it would seem particularly desirable to repeat 
the undertaking on a much larger scale some distance to the south and to 
incorporate the additional groups at first planned, for the biological impli- 
cations of the experiment reach beyond the immediate question of bird 
migration. 

A fuller report of the results, including the histological details of supra- 
renal, thyroid, parathyroid, pituitary and gonad will appear elsewhere. 
It is of sufficient interest to note here, however, that the province-wide 
hunt for crows brought us half a dozen wild birds (all single) that had 
failed to migrate. Histological examination showed gonads of all of these 
to be abnormal, the birds apparently being healthy in all other respects. 

In conclusion, the writer wishes to express his deep gratitude to the 
Directors of the Bache Fund and the Elizabeth Thompson Fund, to 
Colonel John E. Thayer and Dr. John C. Phillips, for generously financing 
the undertaking. He also wishes to express his great appreciation of the 
wide interest taken in the project in Alberta and to thank all and sundry 
who so willingly gave their assistance in one way or another and helped 
to make the undertaking possible. 














